Statistika i osnovna mjerenja 


Uvod 

M. Makek 
2017/2018 


Pravila izvođenja kolegija 


Nastava : 

-2 + 3 sata (predavanja + vježbe) kroz cijeli semestar 

- Laboratorijske vježbe - svaki student radi 3 vježbe 

Kolokviji : 

- 1. kolokvij - uvjet za pohađanje laboratorijskih vježbi 

- 2. i 3. kolokvij - nastavno gradivo 

- Uz određeni broj bodova iz kolokvija moguće je 
oslobađanje od pismenog ispita (pravila na web stranici) 

Uvjet za potpis i pristupanje ispitu : 

- položen prvi kolokvij i pozitivno ocijenjene sve tri 
laboratorijske vježbe 



Prvi kolokvij 


15 . 3 . 2018 . 

u lOh (termin redovite nastave) 
Gradivo : 

- Račun pogreške 

— Grafički prikaz rezultata 

- Metoda najmanjih kvadrata 



Kontakt i informacije 


Sva pravila izvođenja kolegija i dodatne 
informacije dostupne su na web stranici: 

http://www.phy.pmf.unizg.hr/~makek/som/ 

Konzultacije : 

Mihael Makek 
Utorkom 13-14h, soba 326 
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Račun pogreške 

M. Makek 
2017/2018 


Obrada rezultata mjerenja 


Prvi dio kolegija - obrada rezultata mjerenja je uvod 
koji je nužan za pristup laboratorijskim vježbama 

Praktični pristup - primjeri, formule i zadatci 

Sadržaj: 

□Račun pogreške 

□Grafički prikaz rezultata 

□Metoda najmanjih kvadrata 

Cilj : statistički obraditi i pravilno prikazati rezultate 
mjerenja 



RAČUN POGREŠKE 


Vrste pogrešaka pri mjerenjima 
Neovisna mjerenja 
Ovisna mjerenja 

Opća srednja vrijednost i nepouzdanost 
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Mjerenje fizikalnih veličina 


Cilj mjerenja je utvrditi brojčanu vrijednost neke veličine 

Zbog raznih utjecaja rezultat mjerenja x ( -odstupa od prave 
vrijednosti veličine X 

Odstupanje pojedine izmjerene vrijednosti od prave 
vrijednosti se naziva pogreškom mjerenja : Ax=x r X 

Cilj eksperimenta je da izmjerena vrijednost bude što bliže 
pravoj vrijednosti i da pogreška mjerenja bude pravilno 
procijenjena kako bi rezultat bio valjan 



Pogreške mjerenja 


Razlikujemo tri tipa pogrešaka mjerenja 

• Slučajne pogreške 

• Sistematske pogreške 

• Grube pogreške 



Grube pogreške 

Mogu nastati naglim poremećajem u okolini ili 
uređaju - npr. kvar 

Također mogu nastati zbog propusta 
eksperimentatora - npr. krivo očitanje ili zapisivanje 
rezultata 

Mjerenja za koja je ustanovljeno da sadrže grube 
pogreške treba odbaciti iz daljnje obrade ili ponoviti 
u ispravnim uvjetima 



Sistematske pogreške (I) 


Primjer : 

• Mjerimo period titranja matematičkog 
njihala pomoću štoperice 

• Moguće situacije: 

a) Vrijeme na štoperici odstupa od stvarnog vremena za konačan iznos. 
Npr. kad štoperica pokazuje 60 s, u stvarnosti je to 59 s 

-> sva mjerenja vremena su sustavno manja za 59/60 ~ 1.7% 

b) Prilikom pritiska na gumb štoperice potrebno je neko vrijeme rda 
uređaj reagira. 

-> vrijeme reakcije uvijek uzrokuje kašnjenje mjerenja za isti iznos r 

c) Kut očitavanja daje krivi ravnotežni položaj 

-> kut očitavanja uzrokuje pogrešno mjerenje vremena za isti iznos t 








Sistematske pogreške (II) 


Netočnosti mjerenja 


-> npr. pogrešno kalibrirana ili pomaknuta skala 

• Prilikom ponavljanja mjerenja javljaju se u istom smjeru i iznosu 
(reproducibilnost!) 

• Nisu predmetom statističke analize 

• Mogu se ukloniti ili smanjiti na nekoliko načina: 

- poboljšanjem aparature ili tehnike - npr. preciznija kalibracija 

- planiranjem mjerenja — mjerenje se u nekim situacijama može 

organizirati tako da se sistematske pogreške ponište. 

Primjer : kada mjerimo vrijeme štopericom svaki puta kad pritisnemo gumb 
mjerenja kasni za vrijeme r. U slučaju mjerenja perioda njihala isto vrijeme 
kašnjenja je prisutno pri pokretanju i zaustavljanju štoperice, pa se ove dvije 
pogreške poništavaju. 
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Slučajne pogreške (I) 


Primjer : 

• Mjerimo period titranja 
njihala pomoću štoperice 

-> u tom slučaju mjerni 
uređaj čine čovjek + štoperica '"'H+h" 

• Izvršimo 10 mjerenja i dobivamo vrijednosti koje se 
razlikuju 

• Što može dovesti do razlike u rezultatima mjerenja: 

a) Brzine reakcije kod uključivanja/zaustavljanja štoperice 

b) Preciznost štoperice (zaokruživanje decimala) 

c) Nejednako očitavanje ravnotežnog položaja njihala 

nesavršenost mjernog uređaja 

d) Promjene ravnotežnog položaja njihala u vremenu, npr. 
zbog drmanja postolja na kojem se nalazi ili strujanja zraka 

utjecaj okoline 



Mjerenje l 

T [s] 

1 . 

12.5 

2. 

11.9 

3. 

13.0 

4. 

13.1 

5. 

12.6 

6. 

12.2 

7. 

12.7 

8. 

11.6 

9. 

13.4 

10. 

12.8 







Slučajne pogreške (II) 


• Svojstva slučajnih pogrešaka: 

- Ponavljanjem mjerenja dobivaju se različiti rezultati - slučajne 
pogreške različite po iznosu i smjeru 

- Mogu se statistički obraditi (čime ćemo se mi baviti) 

• Uzrok pogreške su nestalni uvjeti mjerenja: 

- Preciznost mjernog uređaja 

- Promjena okoline 

• Smanjenje slučajnih pogrešaka: 

- Usavršavanje mjernog uređaja ili tehnike (npr. automatizacija mjerenja 
u našem primjeru) 

- Izolacija od okoline 

- Ponavljanjem mjerenja i statističkom obradom može se preciznije 
odrediti prava vrijednost fizikalne veličine 

neodređenost rezultata zbog 
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Neovisna mjerenja 


• Izvodimo niz mjerenja neke fizikalne veličine X i 
dobivamo rezultate x p x P x p x n , koji se međusobno 
razlikuju zbog prisustva slučajnih pogrešaka 

• Da bi odredili najvjerojatniju vrijednost mjerene veličine i 
pogrešku mjerenja definiramo sljedeće pojmove: 

> Srednja vrijednost 

> Srednja kvadratna pogreška pojedinog mjerenja 

> Srednja kvadratna pogreška aritmetičke sredine 

> Relativna nepouzdanost 

> Maksimalna apsolutna pogreška 
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Srednja vrijednost 


izračunava se kao aritmetička sredina izmjerenih vrijednosti 


za n mjerenja aritmetička sredina je: 



Uzimamo da je upravo x najvjerojatnija prava vrijednost X 
mjerene fizikalne veličine -> X govori o očekivanoj vrijednosti 
mjerene veličine 
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Srednja kvadratna pogreška mjerenja 

(standardna devijacija) 


• kvadratna odstupanja izmjerenih vrijednosti od srednje 
vrijednosti 



• za dovoljno velik n (~10 mjerenja), m poprima ustaljenu 
vrijednost 

• govori o pouzdanosti pojedinog mjerenja 

• iskazuje rasipanje rezultata kao posljedicu preciznosti uređaja 

mjera preciznosti uređaja 
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Srednja kvadratna pogreška 
aritmetičke sredine (nepouzdanost) 

• Ako izvedemo veći broj mjerenja očekujemo da će fizikalna 
veličina biti preciznije određena 

• Mjera preciznosti rezultata je srednja kvadratna pogreška 
aritmetičke sredine: 

M = 

1 L n 

• M n se smanjuje sa brojem mjerenja proporcionalno ~ 1/Vn 
povećava se preciznost rezultata 

• Govori o pouzdanosti rezultata 


- x > 2 

J n(n - 1 ) 


14 




Pogreške u odnosu na broj mjerenja 
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Maksimalna apsolutna pogreška 

• najveće odstupanje pojedinačnog mjerenja od 
aritmetičke sredine: 

Ax \x ^ilmcuc 

• Ponekad zbog prirode eksperimenta ne možemo 
izračunati M n : 

- ako imamo samo jedno mjerenje 

- ako sva mjerenja daju isti rezultat pa je M n =0 

• Tada procjenjujemo maksimalnu pogrešku Ax 
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Relativna nepouzdanost 

• ako je poznat M n : 

M n 

R = -1- 100 % 
x 

• ako nije poznat M n : 

Ax 

R = —- 100 % 

x 


17 




Rezultat 


• pišemo u obliku: 

x = x + M n 


• ako nije poznat M n : 

x = x ± Ax 
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Zaokruživanje rezultata 


Primjer : 


- pri ponavljanju mjerenja perioda njihala dobili smo srednju 
vrijednost: 7 = 2.453 s 

- pogreška mjerenja iznosi: A7 = 0(2j23s 

- Rezultat pišemo kao: T = 2.5±0.2s 



Prva značajna 
znamenka 


V V 


Pouzdana 


Prva nepouzdana 

znamenka 


znamenka 


Pogreška se u pravilu zaokružuje na prvu znamenku različitu od nule 

Rezulat se zaokružuje na prvu nepouzdanu znamenku - onu koja je 
na istom decimalnom mjestu kao i zaokružena pogreška 
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Primjer 


• Mjerenje perioda T matematičkog njihala 

• Na temelju izmjerenih vrijednosti 
dobivamo: 

T = 12,6 s 
m = 0,6 s 

M 10 = 0,2 s 

• Pa pišemo rezultat mjerenja: 

T = 12,6 ± 0,2 s 


Mjerenje 

T [s] 

1 . 

12,5 

2. 

11,9 

3. 

13,0 

4. 

13,1 

5. 

12,6 

6. 

12,2 

7. 

12,7 

8. 

11,6 

9. 

13,4 

10. 

12,8 
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Broj događaja 


Gaussova raspodjela 


• Ako ponavljamo mjerenje veličine x, dobit ćemo niz rezultata x p x 2 , x 3 ,..., 
x n , koji se međusobno razlikuju zbog prisustva slučajnih pogrešaka 

• Može se pokazati da vrijednosti dobivene ponavljanjem mjerenja slijede 
Gaussovu raspodjelu. (Pri tome zanemarujemo sistematske greške.) 


-3a -2a -a ju a 2a 3a 


1 200 -i 


900 


E00 


300 




V 


II 




99,7% 

-X 


k 


izrnL_ 


Rezultat mjerenja prikazan kao: 
x = x±M y 

♦ I 


srednja vrijednost 


nepouzdanost 


može se opisati pomoću Gaussove 

raspodjele: ( x ~^) 

2ff2 

gdje pretpostavljamo: x = ju 

M =rr/ 
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Ovisna mjerenja 






Tražena veličina Fje funkcija neposredno izmjerenih 
veličina x,, F = f(x 1 , ...,x u ...x^) 

Najvjerojatnija vrijednost je srednja vrijednost: 


F = f (x 1; ..., Xi ,... x n ) 

Ako su veličine x, međusobno neovisne onda je srednja 
kvadratna pogreška veličine F : i- 



• Rezultat pišemo kao: F = F ±M P 
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Ovisna mjerenja 


Primjer: želimo odrediti ubrzanje sile teže mjerenjem perioda 
titranja (7") i duljine niti (/) matematičkog njihala 

• Ubrzanje sile teže je dano relacijom: g = An 2 —7 

i Z 

• Uzmimo da su neovisnim mjerenjima dobiveni rezultati: 

■ / = 0.850 + 0.002 m 

■ T= 1.849 + 0.003 s j 

■ Srednje ubrzanje sile teže je: g = 47 T 2 -=rr 

T L 

■ Srednja kvadratna pogreška je: 



■ Uvrštavanjem dobivamo: g=( 9.82 ±0.04) ms -2 
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Opća srednja vrijednost 

• Izvedeno je n nizova mjerenja iste fizikalne veličine te je za 
svaki niz dobivena srednja vrijednost i kvadratna pogreška: 

= x t ± M t , i = 1... n 

-> primjer fizikalna veličina je određena različitim 
eksperimentalnim metodama 

• Konzistentna mjerenja - ako su razlike za svaki par mjerenja 

| usporedive s bilo kojim M k 

• Nekonzistentna mjerenja - ako su razlike^ — x } \ » M k 
tada zanemarujemo nepouzdanosti M k , a veličine x k 
smatramo nezavisnama te ih tako i analiziramo 
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Opća srednja vrijednost 


Za konzistentna mjerenja definiramo opću aritmetičku 
sredinu: / ” ^ 




Gdje je M nepouzdanost: 

M = 


ZIUMf 2 

Rezultat pišemo u obliku: 

x = x ± M 

Poseban slučaj konzistentnih mjerenja je kada je jedna pouzdanost 
M k znatno manja od svih ostalih. Tada vrijedi m « M k 
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Povećanje preciznosti rezultata 

ili kako smanjiti pogreške? 

• Moguće sistematske pogreške treba reducirati pri planiranju 
eksperimenta - treba razviti i napraviti mjerni uređaj ili 
tehniku koja će omogućiti relativno male sistematke pogreške 
-> bolji uređaj u pravilu znači skuplji uređaj 

• Slučajne pogreške u pravilu se mogu smanjiti ponavljanjem 
mjerenja 

-> dulje mjerenje znači skuplje mjerenje 

• Koliko mjerenja treba napraviti? Toliko da slučajna (statistička) 
pogreška bude manja ili podjenaka sistematskoj 

-> najčešće se pokazuje da je to financijski i vremenski 
najefikasnije rješenje 
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Statističke i sistematske pogreške 


• Primjer mjerenja 

-> Svaka izmjerena točka ima 
pripadajuću statističku 
pogrešku (vertikalne crte) 
i sistematsku pogrešku 
(pravokutnici) 

-> Obje pogreške su 
usporedive 
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Phys. Rev. C 93 014904 
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Grafički prikaz rezultata 

M. Makek 
2016/2017 


GRAFIČKI PRIKAZ REZULTATA 


X-Y graf 
Histogram 


X-Y 


• Pretpostavimo da smo 
mjerenjem fizikalnih veličina 
dobili niz točaka (x^ y ( ) 

• Uobičajeno je da x ona 
veličina koju preciznije 
(neovisno) mjerimo te je 
unosimo na apscisu 

• Izmjerene točke crtamo na X- 
Y graf (eng. scatter plot) 



X-Y graf 


Što uočavamo na grafu: 

a) Linearnost od ishodišta do 
točke A 

b) Nelinearnost od točke A do 
točke C 

c) Rasipanje točaka oko 
zamišljenog pravca u 
linearnom području -> ocjena 
slučajnih pogrešaka 

d) Veliko odstupanje točke B - 
moguće da se radi o gruboj 
pogreški, eventulano se može 
zanemariti. 

Ali: velika odstupanja na 
rubovima grafa poput točke C 
se ne smiju zanemariti jer 
mogu ukazivati na novu pojavu 



X-Y graf 


Prednosti grafičkog nad isključivo numeričkim prikazom: 

- zorno se vidi kako jedna veličina ovisi o drugoj - mogu se uočavati 
zakonitosti, kao npr. linearna ili nelinearna ovisnost na prethodnom 
primjeru 

- vidi se koliko pojedino mjerenje odstupa od trenda - lako se uočavaju 
grube pogreške (poput točke B na prethodnom primjeru) 

Pri tome treba paziti na: 

- Pogodno mjerilo na apscisi i ordinati - graf treba biti popunjen 

- Odabir varijabli za crtanje - ako je moguće linearizirati (o tome kasnije) 

- Označiti fizikalne veličine i mjerne jedinice na svakoj osi i pri tome 
koristiti dovoljno velika slova 

- Kod opisa grafa treba jasno naznačiti što je mjereno i što je prikazano 



X-Y graf 


Interpolacija 

- Može se procijeniti zavisna varijabla y u točki koja nije 
izmjerena, a nalazi se između dvije mjerene točke 

- U pravilu daje ispravne vrijednosti, ali treba uzeti u obzir 
razlučivost 

Ektrapolacija 

- Može se procijeniti zavisna varijabla y u točkama koje se 
protežu izvan mjernog područja 

- U pravilu se izbjegava jer u nepoznatom području može 
doći do odstupanja od trenda - npr. odstupanje od pravca 
između točaka A i C u prethodnom primjeru 



Linearizacija X-Y grafa 


Nakon unošenja točaka u graf uočavamo da ovisnost y(x) nije 
linearna, tj. točke ne leže na pravcu 

Ako očekujemo nelinearnu ovisnost onda pokušamo pravilnim 
odabirom varijabli na x i y osi nacrtati linearan graf. 

Primjer : mjerimo ovisnost prijeđenog puta o vremenu pri jednoliko 
ubrzanom gibanju 

-> očekujemo ovisnost: s(t) = l/2at 2 
-> supstitucija: u=t 2 

-> Crtamo graf: s(u) - u za koji očekujemo da će ovisnost biti 
linearna 



Linearizacija X-Y grafa 


Primjer 1: 

Kvadratna ovisnost s(t) = l/2at 2 linearna ovisnost s(t 2 )=s(u)=l/2a u 
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Linearizacija X-Y grafa 


Ako uočimo nelinearnu ovisnost točaka, a točna funkcijska 
ovisnost y(x) nam nije poznata možemo iskoristiti pravilo 
logaritmiranja: 

- za funkcije tipa y = ax b logaritmiranjem dobivamo: lny = Ina + b\nx 

-> crtamo graf In y- In x u kojem će b odgovarati nagibu pravca a In a 
odsječku na osi y 

- Za funkcije tipa y = ae xb logaritmiranjem dobivamo: ln y = ln a + bx 

crtamo graf In y-x u kojem će b odgovarati nagibu pravca a In a 
odsječku na osi y 


A (nm) 


Linearizacija X-Y grafa 


Primjer 2: 

Nepoznato ovisnost tipa X(L)=AL a logaritmiranje 
linearna ovisnost log X = log A +alog L 


1,6- 
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Histogram 


• Graf koji prikazuje raspodjelu frekvencija 
diskretne ili kontinuirane varijable 

— Diskretne varijable - skup vrijednosti prebrojiv 
— kontinuirane varijable - skup vrijednosti kontinuiran 

• Frekvencija neke određene vrijednosti x, 
varijable X jest broj pojavljivanja te vrijednosti 
u promatranom skupu podataka 

• Relativna frekvencija je definirana kao:/„ =f/N, 
gdje je N ukupan broj podataka 
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Crtanje histograma 


• Za diskretnu varijablu: 

- Odrediti frekvencije ili relativne frekvencije 

- Nacrtati moguće vrijednosti fizikalne veličine na x osi 

- Nacrtati odgovarajuće frekvencije na y-osi 

• Primjer: 
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Crtanje histograma 


• Za kontinuiranu varijablu: 

- Podijeliti interval mjerenja veličine X na konačan broj razreda (mogu 
biti ekvidistantni ili neekvidistantni) 

- Grupirati podatke u razrede i odrediti (relativne) frekvencije 

- Nacrtati razrede na x osi 

- Nacrtati odgovarajuće frekvencije na y-osi 


• Primjer: 


potrošnja goriva 
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Metoda najmanjih kvadrata 


M. Makek 
2016/2017 


METODA NAJMANJIH KVADRATA 


Linearna regresija 
Linearizacija ovisnosti 
Nelinearna regresija 


Linearna regresija 

(veza zavisne i nezavisne varijable linearna) 


Primjer: 

- Mjerimo ovisnost produljenja opruge (/t) o primjenjenoj sili (F) 

— Izmjereni su podaci te je nacrtan graf: 


2.0 n 


H N) 

A(min) 

0.5 

0.37 

1 

0,78 

1,5 

1,08 

2 

1,53 

2,5 

1,89 



* 


« 


F( N) 


i 

3 


Čini se da točke imaju linearnu ovisnost. Kako odrediti pravac 


koji "najbolje" slijedi točke? 





Metoda najmanjih kvadrata 

kod linearne regresije 

Od nekog pravca y=ax+b točka (x i ,y i ) odstupa za s, odnosno 
vrijedi možemo pisati: y,-ox,+ b + s-, 

Princip najmanjih kvadrata: 

Od svih pravaca regresije najvjerojatniji je onaj za koji je 
suma kvadrata odstupanja minimalna 

TI TI 

S(a, b ) = ^ sf = ^ — (ax t + b )] 2 = min 

i= 1 i = i 

dS(a,b) dS(a,b) 

= 0 -— = 0 


odnosno 


da 


db 




Metoda najmanjih kvadrata 

kod linearne regresije 


dS(a,b ) 


da 

dS(a,b ) 


da 



= V 2[y t - (aXi + Z?)](-% £ ) 

i 

= V 2 [~x i y i + (axf + bx t )] 

i 

= 2 ^ + 2a ^ xf + 2£> 









Metoda najmanjih kvadrata 

kod linearne regresije 


dS(a,b) 


db 

dS(a,b ) 


db 



= ^ 2[y £ - (ax t + b)] 

= 2 y £ — 2 a ^ x t — 2b 1 



Vi~ a 



i 

x t — bn — 0 



gi3 h ~ aInXi 

n 







Metoda najmanjih kvadrata 

kod linearne regresije 


Uvrštavanjem izraza (2) u izraz (1) i sređivanjem dobivamo 
izraze za koeficijente pravca: 


n'Zx i y i - (Z*i)(£yi) 

U — 

n 'Zxf - (2 x t ) 2 

II 

3 1 

2 yi-a X Zi ) 






Metoda najmanjih kvadrata 


• Nepouzdanosti koeficijenata slijede iz izraza za a i b te izraza 
za propagaciju pogreške kod zavisnih mjerenja: 



1 

1 

= 

1 

( a > 

■=i > 

2 

a 1 

l 0 -2) 

»5>H 

i= 1 

& 

k i=i ) 

2 ^ 
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Metoda najmanjih kvadrata 


DD 

y, 

ra 

v? 


1 0.5 

0.37 

0.25 

0.137 

0.185 

2 1.0 

0.78 

1.0 

0.608 

0.78 

3 1.5 

1.08 

2.25 

1.166 

1.62 

4 2.0 

1.53 

4.0 

2.341 

3.06 

5 2.5 

1.89 

6.25 

3.572 

4.72 

Z 7.5 

5.65 

13.75 

7.824 

10.365 


0.0 — 
0 


F(N) 


T 

2 


i 

3 
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Metoda najmanjih kvadrata 



a = 0.756 mm/N 

M a = 0.038 mm/N 

L 5 - 

b =- 0.004 mm 

M b = 0.063 mm 

- 


E 10 - 
g 


0.0 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 

0 12 3 

F(N) 
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Linearizacija 


Primjer 1: 

Kvadratna ovisnost s(t) = l/2at 2 linearna ovisnost s(t 2 )=s(u)=l/2a u 

-> dalje se računa s(u) kao da se radi o linearnoj regresiji 



Primjer 2 : linearizacija se vrši logaritmiranjem relacije -> dalje račun za 
linearnu regresiju. 
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Nelinearna regresija 


u slučaju kada zavisna varijabla nelinearno ovisi o nezavisnoj 

primjer-tjerani prigušeni harmonički oscilator: 

A 


y((o) = 


'( 0 »S ~ " 2 )2 + ( 7 )' 


potrebno je odrediti parametre/A, co 0 i z, za koje je suma kvadrata 
odstupanja minimalna: / ' 2 

A 


S(A,it> 0 , t ) = V 


Vi 


+ (f 


)' 


Iz uvjeta minimizacije: 

dS(A, ćo 0j r) 

dA ~~ ° 




dz 


Slijede koeficijenti A, co 0 \ z 









Nelinearna regresija 


• primjer 
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Primjer : matematičko njihalo - mjerenje 
ovisnosti perioda titranja o duljini niti 


https://phet.colorado.edu/sims/pendulum-lab/pendulum-lab en.html 

Zadatak : 

a) Izmjerite period (T) za 5 različitih duljina (L) i napravite 
tablicu L-T 

b) Linearizirajte ovisnost i metodom najmanjih kvadrata 
odredite ubrzanje sile teže 

c) Logaritmirajte ovisnost i odredite a u izrazu T ~ L a 
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Statistika i osnovna mjerenja 


Osnove kombinatorike 

M. Makek 
2016/2017 


OSNOVNI POJMOVI KOMBINATORIKE 


Teorem o uzastopnom prebrojavanju 

Permutacije 

Varijacije 

Kombinacije 


Teorem o uzastopnom prebrojavanju 


Kombinatorijski problem : 

• Želimo popuniti prazna mjesta 1 i 2 

• Imamo n 1 elemenata kojima možemo popuniti mjesto 1 A nl ) 

• Imamo n 2 elemenata kojima možemo popuniti mjesto 2 (B^..., B n2 ) 

• Na koliko se različitih načina mogu popuniti prazna mjesta 1 i 2? 




a 2 

• • • 

A n i 

Bi 

A i B i 

A 2 Bi 

• • • 

Ani^i 

B 2 

AiB 2 

A 2 B 2 

... 

AnA 

• • • 

B„2 

AiB n2 

A 2 B n 2 

• • • 

A n iB n 2 


-> Može se kombinirati svaki A sa svakim B, stoga imamo n 1 n 2 načina 



Teorem o uzastopnom prebrojavanju 


• Pretpostavimo da imamo k praznih mjesta 

• Imamo n 1 elemenata kojima možemo popuniti mjesto 1 A nl ) 

• Imamo n 2 elemenata kojima možemo popuniti mjesto 2 B n2 ) 

• Imamo n k elemenata kojima možemo popuniti mjesto k Y nk ) 

TM: k mjesta možemo popuniti na n 1 n 2 ... n k načina 

-> Dokaz teorema je indukcijom iz prethodnog primjera 


Primjer : 

□ 2 predjela 

□ 3 glavna jela 

□ 3 deserta 


Možemo složiti 18 
različitih jelovnika 




Permutacije bez ponavljanja 


• Pretpostavimo da imamo n različitih elemenata A n koje možemo 
poredati u niz na više načina -> svaki takav niz zovemo permutacija 

• Koliko ima permutacija n elemenata? 

- Ekvivalentno pitanje je: na koliko se načina može popuniti n mjesta s n elemenata? 


□ 1. mjesto možemo popuniti na n načina -> preostaje n-1 elemenata 

□ 2. mjesto možemo popuniti na n-1 načina -> preostaje n-2 elemenata 

□ k. mjesto možemo popuniti na n-(k-l) način 

□ n. mjesto možemo popuniti na n-(n-l) = 1 način 

Prema teoremu o uzastopnom prebrojavanju n mjesta mogu se popuniti na 


= n(n - l)(n - 2)... 3 • 2 ■ 1 ^ ni 


(čitamo: n faktorijela) 


po dogovoru je: 0!=1 



Permutacije bez ponavljanja 


Primjer : koliko riječi možemo složiti iz znakova A, P, N? 

□ Prvo mjesto možemo odabrati na 3 načina: 

■ A_ 

■ P_ 

■ N_ 

□ Drugo mjesto možemo odabrati na 2 načina: 

■ A P _, A N _ 

■ pa_,pn_ 

■ N A _, N P _ 

□ Treće mjesto na jedan način: 

■ A P N , A N P 

■ PAN,PN A 

■ N A P,NPA 

□ Ukupno 3! = 6 načina 



Permutacije s ponavljanjem 


Pretpostavimo da imamo n elemenata A n od kojih je r v r 2 , r k 
identičnih. Koliko ima permutacija takvih n elemenata? 

Primjer : koliko ima permutacija niza AABBB? 


AABBB 

BAABB 

BBAAB 

BBBAA 

ABABB 

BABAB 

BBABA 


ABBAB 

BABBA 



ABBBA 





Za svaki od ovih nizova postoje: 

• 2! identična koji se dobiju 
permutacijama A 

• 3! identična koji se dobiju 
permutacijama B 


n=5, r x =2, r 2 =3 

-> broj različitih permutacija: n! / (rl!r2!) = 5! / (2!3!) = 10 



Permutacije s ponavljanjem 


• Ako imamo r x identičnih elemenata znači da se njihovom zamjenom ne 
dobivaju nove permutacije -- > dobivamo grupe od po r 2 ! identičnih 
permutacija. Broj takvih grupa je nl/rj 

• Od preostalih n-r 1 elemenata imamo r 2 identičnih. Taj skup permutacija 
dijeli se na grupe od po r 2 ! identičnih permutacija. Broj permutacija je 

n!/ r x ! r 2 ! 

• Općenito je onda broj permutacija r 1; ..., r k jednakih elemenata: 


p(n) 


r lt r z ,... t r k = 


n\ 


ri!r 2 !...r fe ! 


• Ovo je poopćenje relacije za broj permutacija 

-> Ako su r lf r 2 ,... r k = 1 (ili 0) onda se ova relacija svodi na slučaj 
permutacija bez ponavljanja 
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Varijacije bez ponavljanja 


Pretpostavimo da imamo n različitih elemenata A n koje želimo 
poredati u niz od r članova, ne dopuštajući ponavljanje jednog te istog 
elementa A k -> svaki takav niz zovemo varijacija n-tog reda i r-tog razreda 

Koliko ima takvih varijacija? 

Zamislimo da treba popuniti niz od r članova s n različitih elemenata (pri 
tome je uvijek r<n): 

- 1. mjesto možemo popuniti na n načina 

- 2. mjesto možemo popuniti na n-1 način 

- r. mjesto možemo popuniti na n-(r-l) način 


Slijedi da je broj varijacija n-tog reda i r-tog razreda: 

(?a (n — r) 

Vs } = n(n — 1)... (n — r -hl) / 



(n — r)! 





Varijacije bez ponavljanja 


• Primjer : 

- Imamo na raspolaganju TI slova 

- Želimo odabrati oznake od 2 slova za registarske tablice 

- Koliko oznaka možemo napraviti? 

-> Obzirom da ovdje razlikujemo poredak slova i ne 
dopuštamo ponavljanje slova, radi se o varijacijama 27. 
reda i 2. razreda 

^V = 271/251=702 


10 



Varijacije s ponavljanjem 


Pretpostavimo da imamo n različitih elemenata A n koje želimo 
poredati u niz od r članova, dopuštajući ponavljanje jednog te istog 
elementa A k -> svaki takav niz zovemo varijacija n-tog reda i r-tog razreda 
s ponavljanjem 

Koliko ima takvih varijacija? 

- 1. broj možemo odabrati na n načina 

- 2. broj možemo odabrati na n načina 

- r. broj možemo odabrati na n načina 

Slijedi da je broj varijacija n-tog reda i r-tog razreda s ponavljanjem: 



• Pri tome r može biti manji, jednak ili veći od n 
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Kombinacije bez ponavljanja 


• Pretpostavimo da imamo n različitih elemenata A n od kojih želimo 
odabrati r članova, ne dopuštajući ponavljanje jednog te istog elementa 

-> Odabrana r-torka se naziva kombinacija n-tog reda i r-tog razreda 

• Za razliku od varijacija ovdje ne pridajemo značaj poretku odabranih 
elemenata -> dvije kombinacije razlikuju se samo ako sadrže različite 
elemente 


Koliko ima takvih kombinacija? 

(li) 

- uzmimo da smo odabrali varijacije V r 

- varijacije grupiramo u skupove koji sadrže iste elemente, ali s različitim 
poretkom -> svaki takav skup će sadržavati r! varijacija 


Slijedi daje broj kombinacija: 


K r — 


(n) 


n! 


r! (n — r)\r\ 



(čitamo: n povrh r) 
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Kombinacije bez ponavljanja 


• Možemo postaviti dva ekvivalentna pitanja: 

- na koliko načina je moguće odabrati r elemenata od njih n? Odg. (n povrh r) 

- na koliko načina je moguće odabrati (n-r) elemenata od njih n? Odg. (n povrh n-r) 


Slijedi: 


(n-r) = (?) 


• Što se lako vidi i matematički: 



n! 

(n-r)! [n — (n —r)]! 



n! 

(n-r)!r! 



13 




Kombinacije bez ponavljanja 


• Primjer loto 7/39: 

- Ukupno 39 brojeva 

- Treba odabrati 7 brojeva 

• ne ponavljaju se 

• poredak nije bitan 

-> radi se o kombinacijama 39. reda i 7. razreda 
-> K = 39!/(7! 32!) = 15 380 937 
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Statistika i osnovna mjerenja 


Teorija vjerojatnosti 

M. Makek 
2016/2017 


Uvod 


• Pokus - bilo koji postupak ili proces koji rezultira opažanjem 

• Ishod - moguć rezultat pokusa (različiti ishodi se međusobno isključuju) 

• Elementarni događaj - pojedini ishod nekog pokusa 

• Prostor elementarnih događaja - skup svih ishoda nekog pokusa 

• Ako je Q prostor elementarnih događaja onda je: 

□ Događaj - svaki podskup od Q 

□ Elementarni događaj - jednočlani podskup od Q 

□ Složeni događaj - višečlani podskup od Q 

□ Sigurni događaj - cijeli Q 

□ Nemoguć događaj - prazan skup 
Primjer : bacanje kocke 

> 6 mogućih ishoda (elementarnih događaja) 

> Složeni događaj - npr. parni ili neparni brojevi 

> Siguran događaj - bilo koji broj od 1 - 6 



Vjerojatnost a priori 

(klasična definicija vjerojatnosti) 


Uzmimo pokus koji završava s konačno mnogo (n) ishoda tj. elementarnih 
događaja, koji su jednako mogući. 


Definicija : vjerojatnost proizvoljnog događaja A je dana omjerom broja 
povoljnih elementarnih događaja n A i ukupnog broja elem. događaja n: 


Iz ove definicije slijedi: 



■ 0 <= P(A) <= 1 (vjerojatnost je između 0 i 1) 

■ n A =0 -> P(A) = 0 (ako niti jedan elementarni događaj ne realizira A, tj. A je 
nemoguć događaj, onda je njegova vjerojatnost jednaka 0) 

■ n A =n -> P(A) = 1 (ako svi elementarni događaji realiziraju A, tj. A je siguran 
događaj, njegova je vjerojatnost 1) 


Nedostaci ovakve definicije su: 

■ pretpostavka o jednako mogućim događajima (tu već pretpostavljamo neku 
vjerojatnost pa je definicija kružna) 

■ Definirana na konačnom skupu događaja 




Vjerojatnost a priori 


• Primjer 1 : bacanje 2 kocke 

- Rezultat svakog bacanja je par brojeva (k^) 

- Ukupno 36 elementarnih događaja čini prostor 
elementarnih događaja 

- Definirajmo složen događaj A kao onaj kod kojeg je suma 
brojeva na obje kocke jednaka 6 

- Elementarni događaji koji realiziraju Asu: 

(1,5), (2,4),(3,3),(4,2),(5,1) n A =6 

P(A) = 6/36 = 1/6 
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Vjerojatnost a priori 


• Primjer 2 : uzorak proizvoda 


— Pretpostavimo da imamo 500 proizvoda 


- Proizvođač garantira da je među njima maksimalno 5% neispravnih 

- Želimo ovo provjeriti uzimajući uzorak od 5 proizvoda 


A. Koja je vjerojatnost da se među tih 5 proizvoda nađu 2 (tj. ~40%) loša? 


1. Na koliko načina možemo odabrati uzorak od 5 proizvoda? ->^5 


( 500 ) 


- (?) 


2. Na koliko načina se može odabrati uzorak u kojem je od 5 proizvoda 2 loša? 

' 25 \ / 475 ' 


*«*»«?»-on — 


Vjerojatnost takvog događaja je: P(A) = 


(V)(Y) 

(?) 


= 0,021 


Koristimo tablicu 
loga rita ma faktorijela 


Obzirom da je vjerojatnost takvog događaja vrlo mala (~2%) moglo bi se zaključiti 
da ako se među 5 nađu 2 loša proizvoda ukupni uzorak mora sadržavati više od 
5% defektnih proizvoda i da garancija proizvođača nije ispunjena. 
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Za prethodni primjer: 

Kako računati s velikim faktorijelama? 


• Može se pojednostavniti: 

a) Korištenjem Stirlingove formule: 



b) Uporabom logaritama faktorijela 
(dano u tablicama) 

c) Logaritmiranjem Stirlingove formule a) + b) 



Vjerojatnost a posteriori 


Primjer : bacanje novčića 

— Ako znamo daje novčić "pošten" onda možemo definirati a 
priori vjerojatnost: 

• Elementarni događaji P i G su jednako vjerojatni p P =p G 

• Slijedi: p P +p G = 1 -> p P =p G =0,5 

— Ako ne znamo da novčić daje jednako vjerojatne ishode 
onda možemo novčić bacati puno puta ( n) i odrediti koliko 
puta se pojavljuje G (n G ), a koliko puta P ( n p ) 

— Tada možemo definirati a posteriori vjerojatnost kao: 


P(G) = 




n 


n 





Vjerojatnost a posteriori 


• Definicija: vjerojatnost događaja A jednaka je relativnoj 

frekvenciji pojavljivanja tog događaja fjA) u nizu od n pokusa: 


fM) = f i 

n 

- vrijedi pod pretpostavkom da su pokusi izvedeni pod 
jednakim uvjetima te da je izveden dovoljan broj pokusa da se 
f r (A) ne mijenja s n, tj. da su relativne frekvencije stabilne 

• Nedostaci ovakve definicije: 

— Što znači dovoljan broj pokusa? 

— Koja je vjerojatnost jednog događaja (ako ne možemo ponavljati pokus 
puno puta)? 
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Suprotna vjerojatnost 


Vjerojatnost da se događaj A ne desi 

- Ne događanje A je događaj koji ćemo označiti sa A 

- Ako j e vjerojatnost za doga đaj A, P(A), onda je vjerojatnost da se A ne 


desi: 


P(A) = 1-P(A) 


- Suprotna vjerojatnost još se označava i sa Q(A) te vrijedi: 


P(A) + Q(A) = 1 


Primjer : Imamo 10 proizvoda od kojih 5 dobrih i 5 neispravnih 

- Ako najednom uzmemo 3 proizvoda koja je vjerojatnost da barem 
jedan od njih bude neispravan - ekvivalentno : koja je vjerojatnost da 
se desi suprotan događaj onom kada su sva tri proizvoda ispravna? 

- Vjerojatnost da imamo 3 ispravna je: P(A) = 3 ) = 1/12 

- Vjerojatnost da nastupi suprotan događaj je: 

P(A) = 1-= 0,917 
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Isključivost događaja 

Definicija : Događaji A i B se međusobno isključuju ako 
istovremeno ne mogu nastupiti oba. 

Analogno : događaji se međusobno isključuju ako se 
istovremeno može pojaviti samo jedan od njih 

-> znači da ne postoji elementaran događaj koji bi 
istovremeno realizirao događaje A i B 

grafički se to može prikazati dis- 
junktnim skupovima točaka u ravnini 

-> skupovi elementarnih događaja koji realiziraju A odn. B 
disjunktni 





Zbrajanje vjerojatnosti 


• Pretpostavimo da se događaji A p A 2 ,A k međusobno isključuju 

• Zanima nas koja je vjerojatnost da se desi ili A ± ili A 2 ,ili A k ? 

• Takav događaj označavamo s A 1 +A 2 +...+A k 

("+" među događajima se čita "ili" i znači da se može desiti samo 1 događaj) 


Vrijedi: 



vjerojatnost da se desi složeni događaj 
jednaka je sumi vjerojatnosti 
komponentnih događaja 


• Dokaz slijedi iz definicije vjerojatnosti (npr. a priori): 





n 


-zs-z™ 
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Zbrajanje vjerojatnosti 


Primjer : bacanje kocke 

- Pretpostavimo da svakom broju na kocki pripada jednaka 
vjerojatnost P=l/6 

- Kolika je vjerojatnost da broj na kocki bude paran? 

->To je vjerojatnost da broj bude ili 2 ili 4 ili 6 

->Vjerojatnost za takav događaj je: 

P(Paran) = 1/6 + 1/6 + 1/6 = 0,5 



Nezavisni događaji 


• Događaji A x i A 2 su nezavisni ako pojavljivanje jednog ne 
mijenja vjerojatnost pojavljivanja drugog 

• Kasnije ćemo pokazati da za nezavisne događaje vrijedi da je 
vjerojatnost istovremenog pojavljivanja događaja A x i A 2 
jednaka umnošku vjerojatnosti pojedinih događaja, tj: 

P(A 1 nA 2 )=P( A 1 )P(A 2 ) 

• Grafički se može pokazati da je skup 
elementarnih događaja koji realiziraju 
\A 1 \A 2 jednak presjeku skupova 
elementarnih događaja za A 2 i A 2 

• Implicitan uvjet je da se događaji A ± i A 2 ne isključuju - drugim riječima 
događaji koji se isključuju ne mogu biti nezavisni 
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Množenje vjerojatnosti 

Pretpostavimo da se događaji A p A 2 ,A k međusobno ne 
isključuju -> mogu se dogoditi istovremeno 

Zanima nas koja je vjerojatnost da se istovremeno dese i A 
A 2 , ■■■, i A k ? 

Takav događaj označavamo s fl A 2 n ... n A k 
("n" se čita "i") 

Može se pokazati da prethodni primjer vrijedi i općenito: 



Vjerojatnost ovakvog složenog događaja je jednaka 
umnošku vjerojatnosti svih komponentnih događaja 







Množenje vjerojatnosti 

Primjer - koincidentno opažanje čestica 

- Uzmimo da želimo opažati reakciju raspada piona: 

7T° -> e + e~y 

- U reakciji se kao produkti pojavljuju elektron, pozitron i foton 

- Efikasnost detektora za opažanje elektrona ili pozitrona je c e+e .= 0,9 

- Efikasnost detektora za opažanje fotona je c y =0,2 

Kolika je vjerojatnost opažanja ovakvog događaja? 

- Efikasnost detektora možemo poistovjetiti sa vjerojatnošću opažanja 

- Ako pretpostavimo da je opažanje svake od čestica nezavisan događaj 
onda je vjerojatnost opažanja ovog događaja jednaka vjerojatnosti za 
istovremeno (koincidentno) opažanje ove tri čestice: 

-> P(n°) = P(e + )P(e~)P(y) =0,162 
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Svojstva nezavisnih događaja 


Tvrdnja: ako su A t i A 2 nezavisni onda su iA 1 iA 2 nezavisni 

• Iz grafičkog prikaza slijedi: 

P(A t n a 2 ) = P(A 1 )~ P(A 1 n a 2 ) 

• Zbog nezavisnosti A t i A 2 možemo pisati: 

p(a 1 n a 2 ) = p(a 1 )- p(a 1 )p(a 2 ) 

= P(A 1 )[1-P(A 2 )] 

= PiAjPfa) 

• Čime je prema definicija nezavisnih događaja dokazana 
tvrdnja 
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Ostali složeni događaji 

• Uzmimo događaje A x i A 2 koji su međusobno ne isključuju 

• Kolika je vjerojatnost da nastupi barem jedan od tih događaja, 
tj. A 1 ili A 2 ili A 1 i A 2 ? 

• Takav složeni događaj označavamo sA 1 \JA 2 


• U specijalnom slučaju kad su A t i A 2 nezavisni događaji, vrijedi: 

P(At U A 2 ) = P(A ,) + P{A 2 ) - P(A 1 )P(A 2 ) 


• Iz grafičkog prikaza slijedi: 

P(A t U a 2 ) = P(A t ) + P(A 2 ) - P(A, n a 2 ) 
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Ostali složeni događaji 

Prethodnu relaciju možemo napisati na drugačiji način: 

— po definiciji: P(^A 1 U A 2 ) = 1 — P{A 1 Uj4 2 ) 

— j 4 i U An znači da se ne dogodi A 2 ili A 2 ili oba 

— Ekvivalentno možemo reći da će se dogoditi događaj suprotan A 
suprotan A 2 : A-± O A 2 

Slijedi: P(A 1 U A 2 ) = 1 - P(A 1 n A 2 ) 

Općenito vrijedi: 

p(A 1 u a 2 u ... uA k ) = i- P(A t nA 2 n ... n A k ) 

A posebno za nezavisne događaje: 

P(A t U A 2 u ... U A k ) — 1 — P(A t )P(A 2 ) ... P(A k ) 






Ostali složeni događaji 


Primjer : 

- A,B,Csu nezavisni događaji koji se međusobno ne 
isključuju 

- U nekom pokusu oni imaju vjerojatnosti: 

P(A)=0.5, P(B)=0.3 i P(C)=0.1 

- Koja je vjerojatnost da nastupi bar jedan od događaja 

P{A U B U C) = 1 - P(a)P(B)p(Č) = 0,685 



Konstantna vjerojatnost 




Događaji čija vjerojatnost se ne mijenja tijekom pokusa 

— Vjerojatnost da događaj nastupi u nekom pokusu: 

P(A 1 ) = P(A 2 ) = - = p 

— Vjerojatnost da događaj ne nastupi u nekom pokusu: 

P(A 1 ) = P(A 2 ) = - = l-p = q 

— Vjerojatnost da događaj nastupi barem jednom u seriji od k pokusa: 

P(A X U A 2 U ... U A k ) - 1 - P(A 1 )P(A 2 ) ... P(A k ) 


(događaji su po definiciji nezavisni jer se događaju u različitim pokusima) 


- Slijedi: P(A 1 U A z U ... U A k ) = 1 - q 


k 
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Konstantna vjerojatnost 


• Primjer : 

- U nekom pokusu događaj A nastupa uz vjerojatnost 
P(A)=p=0,4 

- Koliko pokusa treba učiniti da u njima sa vjerojatnošću od 
95% nastupi barem jedan događaj? 

P = l-q n = 0.95 

-^0,6 n =0.05 

n= 5,86 ~ 6 
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Uvjetna vjerojatnost 


Kolika je vjerojatnost da se dogodi A ako se dogodio B? 


- Označava se sa P(A|B) (čitamo "A ako je B") 

- Skup elementarnih događaja koji realiziraju B onda postaje 
prostorom elementarnih događaja za A| B 


- U tom prostoru AnB elementarnih događaja realizira dog. A 

H “ PM n B) I 

Slijedi da je: P(A\B) = ' , (uvjet P(B)>0) 


Jednako vrijedi: P(B\A) = 


P(A ne) 
P{A), 


Iz čega slijedi: P(A n B) = P(A)P(B\A ) = P(B)P(A\B) 


(vjerojatnost istovremenog zbivanja dva događaja jednaka je produktu 
apsolutne vjerojatnosti prvog i uvjetne vjerojatnosti drugog) 
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Nezavisnost događaja 


• Definicija : događaji A i B su nezavisni onda i samo onda 
kad vrijedi: 

P (A | B) = P(A) 

(Vjerojatnost da se dogodi A ne ovisi o tome da li se dogodio B) 

• Iz relacije za uvjetnu vjerojatnost vrijedi: 

, P(AHB) , N 

p(.m= P(B) - =pw 

• Iz toga slijedi relacija koju smo prije uveli: 

P(AHB) = P(A)P(B ) 
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Potpun sistem događaja 

Neka za događaje A { (\=1,2,...) vrijedi: 

A t * 0, Vi 

Ai D Aj = 0, 

(J Ai = n 

i 

-> Događaji A t čine 
potpun sistem 
događaja 


Vi,;, i * j 










Zakon totalne vjerojatnosti 


Neka /A, čine potpun sistem događaja 


Tada za bilo koji događaj B vrijed 

p(B) = p(Bn^) 


P(B) = Y u P(B\A i )P(A i ) 
















Bayesov teorem 


Neka /A, čine potpun sistem događaja i B neki događaj 
Za bilo koji događaj A k vrijedi: 


P(A k \B) 


PiAnB) 
P(B ) 


Slijedi Bayesov teorem: 


P(A k \B ) 


P(B\A k )P(A k ) 

ZtPiBlAjPiAi) 





Bayesov teorem 


Primjer: kontrola proizvoda 

- Serija proizvoda sadrži 95% ispravnih proizvoda 

- Kontrola proglašava ispravan proizvod dobrim s vjerojatnoću od 98% 

- Kontrola proglašava neispravan proizvod dobrim s vjerojatnošću od 5% 


Koja je vjerojatnost da je proizvod stvarno dobar ako ga je 
kontrola proglasila dobrim? 


- A x = proizvod stvarno dobar, P(A X ) = 0.95 

- A 2 = proizvod stvarno loš, P(A 2 ) = 0.05 

- B = kontrola proglašava proizvod dobrim 


Potpun sustav 
događaja 


- P(AJB)=? 

->P(A 1 \B) = 


P(A 1 )P(B\A 1 ) + P(A 2 )P(B\A 2 ) 


0,997 


-> Nakon kontrole serija će sadržavati 99,7% dobrih proizvoda 




Nedostaci klasičnih definicija 

vjerojatnosti 

• A priori : vjerojatnost proizvoljnog događaja A je dana omjerom broja 
povoljnih elementarnih događaja n A i ukupnog broja elem. događaja n 

• Nedostaci ovakve definicije su: 

■ pretpostavka o jednako mogućim događajima (tu već pretpostavljamo neku 
vjerojatnost pa je definicija kružna) 

■ Definirana na konačnom skupu događaja 

• A posteriori : vjerojatnost događaja A jednaka je relativnoj frekvenciji 
pojavljivanja tog događaja f r (A) u nizu od n pokusa 

• Nedostaci ovakve definicije: 

- Koliki je broj pokusa dovoljan? 

- Koja je vjerojatnost jednog događaja ? 

(ako ne možemo ponavljati pokus puno puta) 
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Aksiomatska izgradnja teorije 

vjerojatnosti 

• Ako poznajemo prostor elementarnih događaja Q za neki 
pokus, svrha definicije vjerojatnosti je da svakom događaju 
A c Q pridruži broj P(A), koji će biti precizna mjera šanse da 
se A ostva ri 

• Definicija vjerojatnosti treba biti općenita i obuhvaćati i 
klasične definicije vjerojatnosti 

• Takvu definiciju temeljenu na aksiomima uveo je Kolmogorov 
1933. godine 

• Objekti u aksiomatskom pristupu su slučajni događaji 
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Kolmogorovi aksiomi 


• Uzmimo da je Q prostor elementarnih događaja 

• Nekom događaju ^jcfi želimo pridružiti vjerojatnost P(A) 

• Pri tome tvrdimo da svako takvo pridruživanje mora 
zadovoljavati (aksiomi): 

1. P(A) > 0, VA 

2. P(il) = 1 

3. a) ako se konačan broj događaja međusobno isključuju (svaki 
par je disjunktan), vrijedi: 


( n \ n 

U ^)= Z PM 

\/ t ■ 


b) ako se prebrojivo beskonačan broj događaja međusobno 
isključuje, vrijedi: 


( go \ on 
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Svojstva vjerojatnosti 


i. 

a) Zbog zatvorenosti skupa elementarnih događaja: A, A G Q 

b) Prema definiciji suprotnog događaja je: A H A = 0 

c) Siguran događaj: P(A U A) = 1 

• Iza), b), c) i trećeg Kolmogorovogaksioma slijedi 

P(A) + P(A) = 1 


Iz prvog aksioma i relacije (I.) slijedi: 


0 < P (A) < 1 


ili« 

Ako u relaciju (I.) umjesto A uvrstimo siguran doga đaj Cl, dobivam 
• Uz postavku da je fi = 0 slijedi: P(0) = 0 





Svojstva vjerojatnosti 


IV. 

— Uzmimo da je A 1 c A 2 
— Onda je A 2 = A 1 U ( A 2 — A 

- Iz trećeg aksioma slijedi daje P(A 2 ) = P{A ± ) + P(A 2 - A ± ) 

— Za bilo koja dva događaja vrijedi i: 

A 1 U A 2 = A 1 U [A 2 — (j 4 1 fl j4 2 )] 

— Pa iz trećeg aksioma i gornjih relacija slijedi: 

P {A U B) = P(A) + P(B) - P {A n B) 

• Sva ova svojstva poznata su i u klasičnoj definiciji vjerojatnosti 
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Podudarnost sa klasičnom definicijom 


• Neka pokus ima velik broj mogućih ishoda - elementarnih 
događaja <&, 

• Tada je vjerojatnost složenog događaja A: P(A) = 

CDjtzA 

• Ako imamo n jednako vjerojatnih ishoda (kao u definiciji a 

priori), onda je vjerojatnost svakog od njih: p(©.) = -,Vt 

n 

• Proizlazi da je vjerojatnost događaja A: n , 

n 

gdje je n A broj elementarnih događaja koji su podskupovi A 
-> to je klasična definicija vjerojatnosti! 
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Statistika i osnovna mjerenja 


Slučajne varijable 

M. Makek 
2016/2017 


Slučajna varijabla 

Q je prostor elementarnih događaja C 0 j (definiranih u 
aksiomatskom, općenitom pristupu) 

slučajna varijabla X je funkcija, koja svakom ishodu pokusa co-, 
iz /^pridružuje neki broj x iz skupa realnih brojeva: 

X\ O -► M 

n.b. slučajna varijabla nije varijabla već funkcija i nije slučajna. 
Naziv potječe od činjenice da je domena od X skup slučajnih 
događaja, a obzirom da se uvriježio mi ćemo ga koristiti. 


Slučajna varijabla 

Ako je Q prostor elementarnih događaja slučajna varijabla 
pridružuje neki realni broj svakom elementarnom događaju 



D je skup svih vrijednosti 
slučajne varijable 













Slučajna varijabla 

Primjer: bacanje 2 kocke 

• Prostor el. događaja Q ima ukupno 36 elemenata co=(iJ) gdje 
i=1..6 i j=1..6 

a) Možemo definirati slučajnu varijablu X koja je jednaka sumi 
brojeva na kockama: X(to)=i+j 

slučajna varijabla X događajima pridružuje vrijednosti 2...12 ili 
skraćeno kažemo slučajna varijabla poprima vrijednosti 2 ... 12 

b) Možemo definirati slučajnu varijablu Y koja je jednaka 
apsolutnoj razlici brojeva na kockama: Y(oo)=|i-j| 

slučajna varijabla Y događajima pridružuje vrijednosti 0...5 ili 
skraćeno kažemo slučajna varijabla poprima vrijednosti 0...5 



Vrste slučajnih varijabli 

Bernoullijeva slučajna varijabla - s.v. koja opisuje događaj koji 
ima samo dva moguća ishoda (Bernoullijev događaj) 

— primjer: bacanje novčića 

Diskretna slučajna varijabla - s.v. za koju je skup mogućih 
vrijednosti konačan ili prebrojiv 

— primjer: bacanje novčića, bacanje kocke, rulet 

Kontinuirana slučajna varijabla - s.v. za koju je skup mogućih 
vrijednosti neki interval na brojevnom pravcu 

— primjer: vrijeme između 2 radioaktivna raspada iz nekog uzorka 



Diskretna slučajna varijabla 

Def. Raspodjelo vjerojatnosti diskretne slučajne varijable 
X definirana je za svaki x E D: 

p(x) = P(X = x) = P({a) E Cl:X(a>) = x}) 

(riječima: vjerojatnost da s.v. poprimi vrijednost x, koja odgovara 
vjerojatnosti da nastupi bilo koji elementarni događaj co kojem je 
pridružena vrijednost x) 


• Pri tome mora biti zadovoljen uvjet da suma svih vjerojatnosti mora biti 
jednaka 1: / p(x) = 1 
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Diskretna slučajna varijabla 

Def. kumulativna funkcija distribucije (funkcija raspodjele) za 
diskretnu s.v. s raspodjelom vjerojatnosti p(xj: 

F(x) = P(X <x) = I p(y)\ 

_ y^x _ 

(riječima: F(x) je vjerojatnost da s.v. poprimi vrijednost manju 
ili jednaku x) 

• pri tome vrijedi: F(— oo) = 0 i F(+oo) = 1 

Za bilo koja dva broja a i b (a<b) vrijedi: P(a<X <b) = F(b)-F(a~) 
gdje je "a-" najveća vrijednost varijable X koja je < a. 



Diskretna slučajna varijabla 

primjeri 


Raspodjele vjerojatnosti 


J 

/ y 

i 


ii 
11 
* 9 
2 9 

2 2 
2 ^ 
2 ž 


I. 

v y 

9 . 9 

y y 

2 9 

9 , 2 

9 2 

Mja 


li 

9 9 

y y 

JM. 


0 1 2 3 4 5 6 7 8 9 10 11 12 13 


I 

M 


] Y//\ 17771 r 


Kumulativne funkcije raspodjele 


0 5 10 15 20 25 30 35 


=*7=-1-1— 


3456789 10 


01 23456789 10 11 12 13 

X 


5 10 15 20 25 30 35 40 

X 


0123456789 10 

X 
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Diskretna slučajna varijabla 

Def. Očekivana vrijednost diskretne slučajne varijable 

• Za diskretnu s.v. X sa skupom mogućih vrijednosti D i 

raspodjelom vjerojatnosti p(x) definirano očekivanje ili srednj 
vrijednost kao: 

E{Xj = ji x = x ■ p(x ) 

xED 



Primjeri : očekivanje diskretne slučajne varijable 


0.2-i 0.2-i 



xEĐ 
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Diskretna slučajna varijabla 


Def. Varijanca diskretne slučajne varijable 

• Za diskretnu s.v. X sa raspodjelom vjerojatnosti p(x) definira se varijanca: 



Varijanca je jednaka sumi očekivanih kvadratnih odstupanja odn. jednaka 
je očekivanju kvadrata odstupanja od srednje vrijednosti 

o x = se naziva standardnom devijacijom ili standardnim 

odstupanjem s.v. X 

Također vrijedi: V(X) = ^ (x 2 — 2xp x + i u 2 )p(x) 


XED 


V(X) = ^ x 2 p(x) — 2j u x ^ xp(x) 4 - j u x ^ p(x) 


xED 


xED 


xED 


V(X) = ^ x 2 p(x) - 2p x p x + 4 |y(jQ = E(X 2 ) - [£(JQ] 


XED 
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Primjer: bacanje 2 kocke 


• Definiramo slučajnu varijablu X kao sumu ishoda bacanja 2 kocke 

• X poprima vrijednosti: x=2...12 

• Raspodjela vjerojatnosti: 



• Očekivanje: E(X)= 7,0 

• Standardna devijacija: cr x = 2,4 
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Funkcija diskretne slučajne varijable 


Očekivanje funkcije g(X) diskretne slučajne varijable X je dano s 

E (gQ 0) = ^ g(x)p(x ) 


XED 


Specijalno za linearnu funkciju vrijedi: 


E(aX + b) = V (ax + b)p(x) = a ^ xp(x) + b ^ p(x) 


XED 


xED 


xED 


— Slijedi: E(aX + b) = aE(X) + b 


- U ovom specijalnom slučaju može se reći da je očekivanje funkcije 
jednako funkciji očekivanja 



Funkcija diskretne slučajne varijable 

Varijanca funkcije g(X) diskretne slučajne varijable X je dana 

V(g(x )) = ^ (#(x) - E(g(X))J p(x) 

XED 

može se pokazati da vrijedi: 

V(gOQ) = E(g 2 (X )) - [E(g(X ))] 2 


Specijalno za linearnu funkciju vrijedi: 

V(aX + b) = a 2 V(X) 

a ctX+b = aa X 


— Odnosno: 




Centralni moment r-tog reda 


Definira se kao: 

M r = £[(X — ju) r ] = ^(x — ju) r p(x) 

_ XED _ 

Iz definicije slijedi: 

□ /w 0 = i 

□ m 3 = 0 

□ M 2 = Y(* - p) 2 p(x) = V(X) 

XED 

-> Centralni moment drugog reda jednak je varijanci 



Pomoćni moment r-tog reda 


• Definira se kao: 



• Iz definicije slijedi: 

□ m 0 =l 

□ m 2 = E(X) = Mx 

pomoćni moment prvog reda jednak je očekivanju slučajne varijable 
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Veza centralnih i pomoćnih momenata 


Može se pokazati da su centralni i pomoćni momenti povezani 
relacijom (dokaz izostavljamo): 

r 



X ^ k C k WM 

k =1 


Npr.: M 2 — m 2 — ml 

-> Što je otprije poznata relacija v(x ) = E(x 2 ) - [fC^)] 2 
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Funkcija izvodnica 

Za diskretnu slučajnu varijablu definiramo funkciju izvodnicu 

g(t ) = ^ e tx p(x ) 

Iz definicije slijedi: 

$'(t = 0) = x ■ p(x) = 

xED 

g"(t = 0) = ^ x 2 ■ p(x ) = m 2 

® Mft = o ) = X* r ' ! ’ w = ” v 

x€D 

pomoću funkcije izvodnice mogu se dobiti svi pomoćni, a 
posredno i centralni momenti slučajne varijable 


Momenti višeg reda 


Služe se detaljnu karakterizaciju raspodjela 
Moment trećeg reda : 

-^koeficijent asimetrije a 3 = — 

(7 

- a 3 = simetrična raspodjela 

- a 3 > 0 nagnuta udesno 

- a 3 < 0 nagnuta ulijevo 

Moment četvrtog reda : 

■^koeficijentspljoštenosti = 


Mi 


o‘ 


- a 4 = 3 normalno spoljoštena raspodjela 

- a 4 > 3 šiljata raspodjela 

- a 4 < 3 široka raspodjela 






Momenti višeg reda 



Sve raspodjele imaju jednako očekivanje i standardnu devijaciju 
-> možemo ih bolje karakterizirati momoću momenata višeg reda 





















Momenti višeg reda 





x 



E(x) = 6 
g= 1,22 



X 


X 
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Statistika i osnovna mjerenja 


Binomna i Poissonova raspodjela 


M. Makek 
2016/2017 


Raspodjele vjerojatnosti 


Želimo matematički i statistički opisati raspodjele varijabli koje 
se javljaju u prirodi u raznim područjima fizike i izvan nje 

♦♦♦ Diskretne raspodjele : Bernoullijeva, Binomna, Poissonova, 
Hipergeometrijska, idr. 


♦♦♦ Kontinuirane raspodjele : Gaussova, Gama, idr. 



Binomna raspodjela 


Primjeri : 

• Bacanje novčića 

— Dva moguća ishoda P i G 
— Za pošten novčić p(P) = p(G) = 0.5 

• Galtonova daska: http://www.geogebra.org/rn/49479 

— Kuglica pada i na svakom koraku može skrenuti lijevo ili desno 
- p(L) = p(D) = 0.5 

Što je zajedničko ovim pokusima? 

-> Imaju dva moguća ishoda 

> Pili G 

> Lili D 

> Uspjeh ili neuspjeh 




Binomna raspodjela 


Binomni pokus: 

1. Sastoji se od n (Bernoullijevih) pokušaja, gdje je n unaprijed zadan 

2. Pokušaji imaju dva moguća ishoda: uspjeh (A) ili neuspjeh (A) 

3. Vjerojatnost uspjeha (p) i neuspjeha ( q=l-p ) je jednaka za sve 
pokušaje 

4. Pokušaji su nezavisni (ishod jednog ne utječe na drugi) 

Binomna slučajna varijabla: 

Za binomni pokus sa n pokušaja definiramo slučajnu varijablu: 

X = broj uspjeha u n pokušaja, X~Bin(n,p) 

X= 0...n 
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Binomna raspodjela 


Def. Binomna raspodjela 

- Za binomni eksperiment i binomnu slučajnu varijablu X~Bin(n,p), 
njenu raspodjelu vjerojatnosti nazivamo binomnom raspodjelom i 
označavamo: b(x; n,p) 

Uzmimo da niz od n Bernoullijevih eksperimenata završi: 

AAAAAA . AA 

-A se pojavljuje x puta, A n-x puta 

- Koja je vjerojatnost jednog takvog ishoda? -> p x q n ' x 


n 


- Koliko ima takvih istih ishoda (koji su svi jednako vjerojatni)? -> 

- Vjerojatnost da se pojavi bilo koji ishod koji ima x uspjeha je: 



za x=0,l,2,...n 
b(x; n, p) = 0 za sve ostale 


-> vjerojatnost da u seriji od n pokusa događaj A nastupi x puta 




Binomna raspodjela 


• Primjer : Galtonova daska 


b(x; 

b(x; 



slučajna 

varijabla 

X 

putevi 

vjerojatnost 

0 

LLLL 

q 4 

1 

LLLD 

LLDL 

LDLL 

DLLL 

pq 3 

2 

LLDD 

LDLD 

LDDL 

DLLD 

DLDL 

DDLL 

p 2 q 2 

3 

LDDD 

DLDD 

DDLD 

DDDL 

p 3 q 

4 

DDDD 

p 4 












Binomna raspodjela 


Rekurzivna formula 

- Pomaže u računu, posebice kod velikih faktorijela 

- Koja je vjerojatnost da varijabla poprimi vrijednost x-l? 

b(pc — 1 ;n, p) = ^ ^p x ~ 1 q n ~ x+1 

— Iz omjera b(x)/b(x-l) dobivamo rekurzivnu formulu: 

n — x + 1 p 

b(pc; n, p) = - b(x — 1; n, p) 

X 


Vjerojatnost za bilo koji x možemo dobiti rekurzijom 

b(0; n,p) = q n 




Binomna raspodjela 


Primjer korištenja rekurzije : za koju vrijednost x je b(x) 
maksimalan? 

• Koristimo rekurzivnu formulu 
Mora biti zadovoljeno: 

b(x — 1; n, p) < b(x ; n,p) > b(x + 1; n,p) 


Nakon raspisivanja (domaća zadaća) po rekurzivnoj formuli 
slijedi: 


np-q< x max <np + p 
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Binomna raspodjela 


Vrijedi binomni poučak: 


n 


(a + b) n = ^ a k b n 


-k 


k=Q 


Funkcija raspodjele diskretne slučajne varijable sa binomnom 
raspodjelom: 


X 


F(x) = 


-\I0 


p k q n k '0 < x < n 


k =o 


1 , 


x = n 
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Binomna raspodjela 


Očekivanje diskretne s lučajne var ijable sa binomnom 


raspodjelom je: E(X) = np 


x=0 


n 


a :=0 


To možemo pokazati na sljedeći način: 

- Definiramo funkciju g(t) = (q + pt) n 

- Razvoj po binomnom poučku: n 

(q 4- pt) n = ^ p x q n ~ x t x = ^ b(x) 

- Deriviranjem po t dobivamo: x=0 n 

np(q + pt) n_1 = ^ x b(x)t x ~ 1 

- Za t=l: x n° 

np(q + p) n_1 = ^ x b(x) 

- Što po definiciji daje: 

np = E(X) 



Funkcija izvodnica 


Općenito za diskretnu slučajnu varijablu definiramo funkciju 
izvodnicu: 


g(t) = ^ e tx p(x) 


XED 


Iz definicije slijedi: 

' r (t = 0) = ^ x ■ p(x) = m 1 


8 


xED 


8 


"(t = 0) = ^ x 2 ■ p(x ) = m 2 


XED 


» W(t = 0 ) = X* r ' 


x€D 


pomoću funkcije izvodnice mogu se dobiti svi pomoćni, a 
posredno i centralni momenti slučajne varijable 


Binomna raspodjela 


• Drugi način da odredimo očekivanje (a i varijancu) je pomoću 
funkcije izvodnice 

• Funkcija izvodnica za binornnu raspodjelu: 

g x (t) = E(e tx ) = ^ e tx p x q n ~ x 

x = Q 
n 

= V Q ( pe t ) x q n ~ x 

x=0 

• Prema binomnom poučku slijedi: 


gx(t) = (pe f + q) n 
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Binomna raspodjela 


Prva derivacija funkcije izvodnice: 

g' x (t) = npe t (pe t + q) n_1 
g' (0) = npe°(pe° + q) n ~ 1 = np 



E(X ) = np 



Druga derivacija funkcije izvodnice: 

g'x(t) = npe t (pe t + q ) n ~ 1 

+ n(n — 1 )(pe t ) 2 (pe t + q) n ~ 2 

g'x( 0) = np + n(n — l)p 2 (= m 2 ) 





Binomna raspodjela 


• Varijanca je jednaka centralnom momentu drugog reda 

• Otprije imamo: M 2 = m 2 - ml 

• Slijedi Varijanca diskretne slučajne varijable sa binomnom 
raspodjelom: 

V(X ) = np + n(n — 1 )p 2 — (np) 2 
= np( 1 + (n — l)p - np) 

= np(l — p) 


7(Z) = npg 


Treba još primijetit da se za n=l binomna raspodjela svodi na 
Bernoullijevu, što je za očekivati buduć da ona opisuje 
Bernoullijev pokus koji se ponavlja n puta! 
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Binomna raspodjela 


Primjer : kontrola uzoraka 

• Neki uređaj izrađuje 8% neispravnih proizvoda. Proizvodi se 
pakiraju u kutije od 50 komada. 

Koja je vjerojatnost da se u kutiji nađe 2 defektna komada? 

- n=50, p=0.08 



Koji je najvjerojatniji broj defektnih komada u kutiji? 

- n=50, p=0.08 

- np-q <= x M <= np+p -> 3.08 < = x M <= 4.08 
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Poissonova raspodjela 


• Uzmimo da je vjerojatnost uspjeha pojedinog Bemoullijevog 
pokušaja u binomnom pokusu jako mala, te taj pokus 
ponavljamo jako puno puta (n », p « ) 

pri tome očekivana vrijednost konstantna E(x)=np= const. 


• Primjer : raspad radioaktivnog uzorka: 

— Vjerojatnost pojedinog raspada je iznimno mala 
— Uzorak sadrži jako puno izotopa koji se nezavisno mogu raspasti 
-> Srednji broj raspada u nekom vremenskom intervalu je konstantan. 
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Poissonova raspodjela 


Kao granični slučaj binomne: 

b(x; n,p) = ( ) p x q n ' 


1—X 


n(n — 1)... (n — x 4- 1) 

v I 

J\f ■ 

n(n — l)(n — x + 1) X 


y 1 

1\ 2 


n 

x — 1 


ri x n n-x 
r 4 


A' 

n. 


n-x 


x 




x\ 


r , A* / A\ n 

p(x) = lim o(x) = — lim 1-1 

n->oo x! n->oo \ 72/ 


n 



X 

n 


) 


X 



X x . 

p(x) 

- 

r 1 

■ 
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p(x:a) 


Poissonova raspodjela 


• P.r. je potpuno određena parametrom X 





II 

O 

Oi 




■ A =1 

0.6 - 


A 

2=2 




■ 2=3 

- 



II 

CJi 

0.4 - 





0.2 


0.0 





X 
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Poissonova raspodjela 


Očekivanje je prema definiciji: 

00 


E(X) = £ 


x 


X 


x — e 

v! 

/V ■ 

x=0 00 


-X 


= Xe~ A 


E(X) = X 


I 

x=0 


X 


X —1 


(x — 1)! 


Prema definiciji varijance također izlazi da je: 


V(X) = X odn. 



Normiranost raspodjele slijedi iz definicije: 


QO QO 

X x 


e = 


x=0 


x=Q 


Xe A e A 
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Poissonova raspodjela 


• Funkcija izvodnica: 

GO 

Z X x 
e tx e -A_ 

x=0 


e -AgAe c _ g A(e t -l) 


• Pomoćni momenti: 

g'(t ) = Ae t e^ et ~ 1 ^ m 1 = A 

g"(t ) — Ae c ( 1 + /le t )e /, '( et-1 ) m 2 — A 2 + A 


E(X) i V(X) možemo elegantno dobiti preko funkcije izvodnice 
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Poissonova raspodjela 


Rekurzija 

- Iz omjera p(x)/p(x-l) slijedi: 

I 

p(x) = — p(x — 1) 

X 


Najvjerojatnija vrijednost 

- Iz uvjeta: pO M “ 1) ^ P(*m) ^ P(*m + 1) 

- Slijedi: U. - 1 < x M < A 




Raspodjela dvije nezavisne 
Poissonove varijable 

Za dvije slučajne varijable X i Y s Poissonovim raspodjelama 
p(x,m) i p(y,n) vrijedi da je raspodjela slučajne varijable Z=X+Y 
jednaka p(z, m+n) 

Dokaz: p(z = z) = Y P(X = k)P(Y = z-k ) 


Zj 

=z 


k =0 

n 

- e -n- 


z—k 


k\ 


k=0 

0 -(m+n) 


z\ 

e -(m+ri) 


10 


(z — k)l 


m k n z k 


k =o 


z\ 


(m + n) z — p(z, m + n) 
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Poissonova raspodjela 

(primjene) 


Primjer : Poissonova raspodjela kao granični slučaj binomne 

• Serija proizvoda sadrži 4% loših 


• Kolika je vjerojatnost da u seriji od 10 proizvoda nadjemo 


najviše jedan loš proizvod? 


- N=10 

- P=0.04 


b(0;10,0.04) + 
b(l;10,0.04) = 
0.9418 


- Np = 0.4 


P(0;0.4) + 
' p(l;0.4)= 
0.9384 


Poissonova r. je dobra 
aproksimacija binomne! 



x 
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Poissonova raspodjela 

(primjene) 


Primjer : Poissonova raspodjela za proučavanje rijetkih događaja 

• Konstanta raspada nekog izotopa je X=7,3 10 10 s 1 

• Uzorak sadrži 4,8 10 22 radioaktivnih atoma 

• Treba odrediti raspodjelu vjerojatnosti k raspada u 1 ms. 

• Zakon radioaktivnog raspada daje broj atoma u vremenu: 

N(t) = N 0 e~’ ,T 






I aktivnost uzorka: 


a = 


dN 


dt 


N 


Slijedi da je vjerojatnost k raspada u intevalu At: 


P k (At) = b(k; N, — —) = p(t, aAt ) 

t N 

Vjerojatnost događaja proporcionalna je veličini intervala 
-> Takve procese nazivamo Poissonovim procesima! 
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Prilagodba empiričkim podacima 

• Primjer Poissonova raspodjela 

• Mjeren je broj beta raspada x iz nekog uzorka u intervalu od 10 sekundi. 
Broj slučajeva u kojima je emitirano točno x čestica je označen kao/ x 

-> smatramo da bi broj opaženih čestica trebao slijediti Poiss. raspodjelu. 


X 

t 

P(x) 

t. 

0 

57 

0,0209 

54,4 

i 

203 

0,0807 

210,5 

2 

383 

0,1562 

407,4 

3 

525 

0,2015 

525,5 

4 

532 

0,1949 

508,4 

5 

408 

0,1509 

393,5 

6 

273 

0,0973 

253,8 

7 

139 

0,0538 

140,3 

8 

45 

0,0260 

67,9 

9 

27 

0,0112 

29,2 

10 

10 

0,0043 

11,3 

11 

4 

0,0015 

4,0 

12 

2 

0,0005 

1,3 


>Očekivanje raspodjele možemo aproksimirati 
aritmetičkom sredinom uzorka: X = x =3,87 
>Teorijske frekvencije računamo prema: ftx = Np(x;X) 



X 






























Druge raspodjele diskretnih varijabli 


• Hipergeometrijska raspodjela 

- Populacija od N elemenata (konačna populacija) 

- Svaki element može biti uspješan (A) ili neuspješan (A) 

- Postoji M uspješnih elemenata 

- Izvlačenje bilo kojeg uzorka od n elemenata je jednako vjerojatno 

- Slučajna varijabla X=brojuspješnih elemenata u uzorku 

- h(x)=broj mogućih uzoraka s x uspjeha/ukupni broj mogućih uzoraka 


x = 0,1,2,..., min (n, M) 

Ti/ 


/i(x;n,M, JV) = 


("X 


N—M 

n—x 


/w\ 
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Druge raspodjele diskretnih varijabli 


Hipergeometrijska raspodjela 

- Očekivanje E(X) = n — 


- Varijanca V(X) = n—— _——_- 

N N N -1 

- Uz definicije: p=M/N i q=(N-M)/N slijedi: 

E(X) = np V(X) = npq 

N ~ 1 

Primjer : skup ima 50 proizvoda od kojih 40 dobrih i 10 loših. Slučajno 
odaberemo 5 proizvoda. Koja je vjerojatnost da se u uzorku nađe x dobrih 
proizvoda? -> h(x; 5,40,50) 



Druge raspodjele diskretnih varijabli 


• Pascalova raspodjela 

- Niz nezavnisih pokušaja 

- Svaki pokušaj je Bernoullijev pokus koji može biti uspješan (A) ili 
neuspješan (A) 

- Vjerojatnost uspjeha p jednaka je za sve pokušaje 

- Pokušaji se izvode sve dok se ne opazi r uspjeha (r unaprijed određen) 

- Slučajna varijabla X = broj neuspjeha koji prethode r-tom uspjehu 

- Za razliku od binomne broj pokušaja je varijabilan, a broj uspjeha je 


zadan 


nb(x; r,p) = P(r — lAux+r— 1 po kušaj d)P (uspjeh) 
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Čebiševljev teorem 


Za slučajnu varijablu X sa raspodjelom vjerojatnosti P(x) i očekivanjem 
E(X)=ju i varijancom V(X)=d* vrijedi teorem: 

Vjerojatnost da X poprimi vrijednost izvan intervala (p-e, p,+8)je manja ili 
jednaka a 2 /e 2 , za bilo koji e>0. 



|Ll-8 p p+S x 

Dokaz : po definiciji je varijanca jednaka: a 2 = y](x i - ^yp(x i ) 
Posebno je za sve x iz osjenčanog intervala: 1 

a 2 > ( x ; “ v) 2 P(xd 

Obzirom da je: \x t — // 1 > 6 onda I 

vrijedi: 

< 7 2 > 6 2 

\Xi-fl\>€ 


^ P(xd 


► 


Budući da se X: međusobno isključuju 


vrijedi: i P(Xi) = P{\x-n\ > €} 
|^,-/i|>fc Odavde slijedi tvrdnja 
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Zakon velikih brojeva 


• Za Bernoullijev događaj s vjerojatnošću P(A)=p, koji se ponavlja u n 
pokusa, vjerojatnost slijedi binomnu raspodjelu sa očekivanjem ju=np i 
varijancom o 2 =npq. Empirijski je relativna frekvencija tog događaja: / 7 
(gdje je x broj nastupa u n pokušaja) 

Bemoullijev teorem: ako broj pokušaja teži u beskonačno onda relativna 
frekvencija teži po vjerojatnosti stvarnoj vjerojatnosti događaja: 

(za bilo koji c) 


TX 


Dokaz 


lim P{\f rx - p\ > e} = 0 


n^co 


p{\frx - pl > e} = P fl- - p > e] = P{|* - np\ > ne] 

Un J <7 2 


Prema Čeb. Tm (uz np=|Li): P{\x — np\ > ne] < 


n 2 e 2 


H\f rx -P\ ž <0 < 


o 


npq 


n 2 e 2 n 2 e 2 

Za lim(n-> beskonačno:) P{\f rx ~ P\ ^ 6 ) = 0 


P{\frx~P\ ^ e} < 


M 

ne 2 


(Q.E.D.) 


x 

n 
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Statistika i osnovna mjerenja 


Kontinuirane raspodjele 

M. Makek 
2016/2017 


Kontinuirana slučajna varijabla 

Definicija : Slučajna varijabla X je kontinuirana ako skup njenih 
vrijednosti čini interval brojeva, tj. ako je za A i B ( A<B) svaki 
x € [A,B] moguć. 


• Kontinuirana slučajna varijabla će biti reprezentacija 
opservabli poput mase, duljine, vremena, itd. 

• Većina karakterističnih veličina (očekivanje, varijanca, idr.) 
kontinuirane slučajne varijable može se dobiti iz veličina 
diskretne slučajne varijable uz prelazak sa sume na integral 



Kontinuirana slučajna varijabla 

• Primjer: mjerimo dubinu jezera na ~10000 mjesta 


zaokruženo 
na 1 m 



X 



X 


zaokruženo 
na 1 cm 


zaokruženo 
na 1 dm 



Prijelaz s diskretne 
na kontinuiranu: 


p(x)^f(x)-dx 



kontinuirana 

raspodjela 


co 

S -J 


x<eD 


—oo 
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Funkcija gustoće vjerojatnosti 

Definicija : funkcija gustoće vjerojatnosti kontinuirane slučajne 
varijable X je funkcija f(x )takva da za bilo koja dva broja a<b 
vrijedi: 


Iz definicije slijedi da je: 

P(X=x 0 ) = 0, za svaki x 0 ! 

Vjerojatnost da X poprimi vrijednost u intervalu [a,b] dana je 
površinom ispod funkcije gustoće vjerojatnosti u tom intervalu 

• Pri tom ef(x) zadovoljava svojstva: 

1 . f(x ) > 0, Vx 

CQ 

2. f f{pc)dx = 1 



— co 


4 




Funkcija raspodjele 


Funkcija raspodjele (kumulativna funkcija distribucije) diskretne 
slučajne varijable X definirana je: 



• Općenito vrijedi: 


P(a < X < b) = F(b ) - F(a)-> distribucijska 


funkcija ima važnu ulogu kod kontuiruiranih si. var. jer pomoću 
nje jednostavno dobivamo vjerojatnost 




Iz distribucijske funkcije jednostavno dobivamo funkciju 
gustoće vjerojatnosti iz: 



dF(x) 


dx 





Funkcija gustoće vjerojatnosti 

Primjer : uniformna ili pravokutna raspodjela u intervalu [A,B] 

( 1 


f(x) = \ B -A 

l o 


A < x < B 
inače 







Funkcija raspodjele 


• Primjer : uniformna ili pravokutna raspodjela u intervalu [A,B] 


F(x) = 


0 , 

x — A 

b^Ta' 

1 , 



X 
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Očekivanje i varijanca kontinuirane 

slučajne varijable 

Očekivanje kontinuirane slučajne varijable X definira se: 



Varijanca kontinuirane slučajne varijable X je: 


V(X ) = crj = I (% — /J.) 2 f(x)dx = E\(X — fi) 2 ] 

— co 

Može se pokazati da i ovdje vrijedi: 

V(X) = E(X 2 ) - [E(X)] 2 
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Očekivanje i varijanca kontinuirane 

slučajne varijable 


Primjer : uniformna ili pravokutna raspodjela u intervalu [A,B] 


1 


f(x ) 


B — A 
0 


A < x < B 
inače 


A + B 

Očekivanje: E(X) = —-— 


Varijanca: 


V(X) = 


(A - B) 2 
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Funkcija izvodnica i momenti 


Definicije : 


funkcija izvodnica: 

co 

g(t) = j 

— OO 

pomoćni moment: 

oo 

f 


7?V = 1 X 

— co 

centralni moment: 

co 


e tx f(x)dx 


M r = J(x- iiYf(x)dx 


— oo 
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Očekivanje i varijanca funkcije 
kontinuirane slučajne varijable 




Za funkciju h(X) kontinuirane slučajne varijable X s gustoćom 
vjerojatnosti f(x), očekivanje je: 


E(h(X)) = fi h ( X ) = j h(x)f(x)dx 


• Za funkciju h(X) kontinuirane slučajne varijable X s gustoćom 
vjerojatnosti f(x), varijanca je: 

V(h(X)) = E[h(X) 2 ] - [E(h(X))f 


Specijalno : linearna funkcija h(X)=aX+b 

— Očekivanje: E(aX+b)=a E(X)+b 

- Va rija nca : V(aX+b)=a 2 V(X) 
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Normalna (Gaussova) raspodjela 


Slučajna varijabla X ima normalnu 
raspodjelu ako je njena funkcija 
gustoće vjerojatnosti: 



gdje x može biti iz (— 00 , 0 °) 

Slučajnu varijablu koja ima 
normalnu raspodjelu 
označavamo sa: 

X~N(pi,cč) 


c = 0.5 



X 
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Normalna (Gaussova) raspodjela 


Normiranost raspodjele slijedi iz definicije: 


e 2 <t 2 dx 

— oo 


e 2 ° 2 clu = 2 

— CQ 


Slijedi: 





uz supstituciju: u = x — fi 
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Normalna (Gaussova) raspodjela 


Očekivanje - prema definiciji: 


co 


co 


E(X) = j x f{x)dx = j 


x 


— co 


— co 


i 


ay/2n 


2c 2 dx 


• Uz supstituciju: 

• Slijedi: 


u = 


x — pi 


a 


E{X) = /i 


Varijanca- prema definiciji se može pokazati da slijedi: 

V(X) = a 2 
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Normalna (Gaussova) raspodjela 


Standardizirana raspodjela 

• Uvodimo substituciju: 


u = 


x — n 


a 


Funkcija gustoće vjerojatnosti standardne normalne 
raspodjele: 



-> Normalna raspodjela s očekivanjem 0 i varijancom 1 U~N(0,1) 

• Veza s "običnom" Gaussovom raspodjelom: 

1 

f(x) = - (p (u) 
o 

• Kumulativna funkcija standardizirane normalne raspodjele: 

ll 

= J ~J^ e 2 Vrijednosti dane u tablicama! 


— co 
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Normalna (Gaussova) raspodjela 


Vjerojatnost da x poprimi vrijednost u intervalu (x 1; x 2 ) je po 


■ r* • • • • 

definiciji: 


x 2 


l r 

P(x-t < X < X?) = ,_ 

\[2na J 


(x-d) 2 

e 2 <j 2 dx 


*i.. 


Prelaskom na standardiziranu raspodjelu, tj. uvođenjem 

_—— vjerojatnost postaje: 

a u 2 


supstitucije u = 


1 r _h 

P(x -, < x < x 2 ) = —== e 2 

\f2n J 


V2tt 

Koristimo distribucijsku funkciju: 


2 

2 ~du 


U 


P(x 1 < x < x 2 ) = <t>(u 2 ) — i>(ui) 
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Normalna (Gaussova) raspodjela 


Vjerojatnost - primjer 

• Uzmimo normalnu raspodjelu p=50 i g=2 

• Kolika je vjerojatnost da X poprimi vrijednost između 46 i 53 

P(46 < X < 53) =? 

x — fi 

• Uz substituciju u = -: 

a 

P( 46 < X < 53) = P(-2 < U < 1,5) 

= 0(1,5) - 0(—2) 

• Tablice: http://www.stat.ufl.edu/~athienit/Tables/Ztable.pdf 

P = 0,9332-0,0228=0,9104 
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Normalna (Gaussova) raspodjela 


• Vjerojatnost da X 
poprimi vrijednost: 


Interval 

P(x) 

ji+a) 

0,683 

(ju-2a, ju+2a) 

0,954 

(ju-3a, ju+3a) 

0,997 


• Simetričnost a 3 =0 

• Spljoštenost a 4 =3 
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Značaj normalne raspodjele 

1. Objašnjava najveći broj statističih opažanja 

2. Opisuje raspodjelu slučajne pogreške mjerenja 

3. Općenito opisuje pojave gdje postoji velik broj slučajnih 
utjecaja -> pokazuje sa da je ukupna raspodjela veličine 
podvrgnute tim utjecajima, Gaussova raspodjela 
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Aproksimacija binomne raspodjele 

Gaussovom 


• Pokazuje se da kod binomne raspodjele: 

— koeficijent asimetrije a 3 -> 0 kad n-> oo « 3 - <t r — 

yjnpq i_6 pq 

- koeficijent spljoštenosti a 4 -> 3 kad n-> oo « 4 =-+ 3 

H npq 

-> to su upravo svojstva normalne raspodjele 
Binomna teži u normalnu kad n oo 

• U praksi se uglavnom pokazuje da je aproksimacija 
Gaussovom raspodjelom dobra ako su zadovoljeni uvjeti: 


l. 


1 n 

-< p < - 

n+1 n +1 


2. npq > 9 
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Aproksimacija binomne raspodjele 

Gaussovom 


0,193 


Primjer : koja je vjerojatnost da u 12 bacanja novčića ispadne 
najmanje 4 a najviše 7 glava? 

• Binomna raspodjela: 

□ n=12 

□ p=l/2 

->P (4 < X < 7) =0,733 

• Aproksimacija Gaussovom: 

© > 1/13 < 1/2 < 12/13 
\ > npq=3<9 

ju = 6, a 2 = 3 

P(3.5 < X < 7.5) =0,733 


0.25 


0.20 - 


0 . 15 - 


0 . 10 - 


0 . 05 - 


0.00 



X 
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Aproksimacija Poissonove raspodjele 

Gaussovom 


• Ako X ima Poissonovu raspodjelu X~P(A.) sa velikim X (>>20) 
može se dobro aproksimirati Gaussovom X~N(>,,A,) 


Primjer : Geigerov brojač mjeri raspa¬ 
de sa srednjom brzinom od 25 s 1 

Koja je vjerojatnost da opazimo naj¬ 
manje 23 i najviše 27 raspada/s ? 

PoissonX~P(25): 

P(23 < X < 27) = 0,3826 

Gauss X~(25,25): 

P(22.5 < X < 27.5) = 0,3830 
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Obitelj gama raspodjela 


Definiramo T-funkciju: 

• Svojstva : 

1. Va > 1 -> r(a) = (a — l)r(a — 1) (rekurzija) 

2. Vn £ M -> T(n) = (n — 1)! (poopćenje faktorijela!) 


CQ 

r(a) — f x a ~ 1 e~ x dx, 

a > 0 

J 

0 









Standardna T-raspodjela 


Xje: 


1 

H 

1 

i 

* 

x > 0 

f (x; a ) = j 

I r (a) ' 

1 

10 , 

X < 0 


raspodjele a>0 
















Standardna T-raspodjela 


Svojstva - osnovni uvjeti raspodjele vjerojatnosti: 


1. f(x; a ) > 0, 


oo 



— co 


Vx 


- [ x a 1 e X dx 

r (a)J 


r (g) 
r(a) 



Očekivanje: E('Xj=a 


Varijanca: \/(X)=a 




Općenita T-raspodjela 

• Funkcija gustoće vjerojatnosti: 

p«r(a) 

0 , 

gdje su a, p >0 

• Očekivanje: E(X)=afi 

• Varijanca: V(X)=a/3 2 
Za P=1 općenita standardna T-raspodjela 


x > 0 
x < 0 
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Statistika i osnovna mjerenja 


Dvodimenzionalne raspodjele i korelacija 


M. Makek 
2016/2017 


Višedimenzionalne raspodjele 

Primjer : u eksperimentu sudaramo 2 protona 

- Prije sudara oni imaju energiju i impuls E v p v E 2 , p 2 

- Nakon sudara nastaje mnoštvo različitih čestica s energijama i 
impulsima E k ',p k ' 

- Da bismo odredili energiju i impuls čestice u konačnom stanju mjerimo 
npr. njen impuls u magnetskom spektrometru i energiju u kalorimetru. 

- U tom slučaju energija i impuls su kontinuirane slučajne varijable 

- Pitanja: da li su nezavisne? Ako ne, kako su one povezane? Kako su 
povezane sa energijama i impulsima ostalih čestica? 

• U eksperimentima se često opaža veći broj varijabli, koje 
mogu biti neovisne ili povezane, o čemu će ovisiti analiza 
rezultata i račun pogreške 



Dvodimenzionalne raspodjele 


Uzmimo dvije slučajne varijable X i Y definirane na istom 
prostoru elementranih događaja 

Def : združena raspodjela vjerojatnosti p(x,y) za diskretne 
slučajne varijable X i Y je vjerojatnost da istodobno X poprimi 
vrijednost x i Y poprimi vrijednost y: 

P(x,y) = P(X=x i Y=y) 

Def : Rubne raspodjele vjerojatnosti za varijable X i Y 
označavamo s p x (x) i p Y (y), a dane su izrazima: 

Px(x) = V p(x,y ) Py(y) = 2_, V(x,y) 

yEDy x ^ D X 

vjerojatnost da pojedina varijabla iz združene raspodjele 
poprimi neku vrijednost bez obzira na drugu varijablu. 




Dvodimenzionalne raspodjele 


Primjer: studenti prve 

Ocjena 

Muško 

Žensko 

UKUPNO 

godine razvrstani po 

1 

10 

2 

12 

ocjeni iz statistike 

2 

13 

3 

16 

(varijabla X) i po 

3 

15 

7 

22 

spolu (varijabla Y) 

4 

18 

7 

25 


5 

4 

3 

7 


UKUPNO 

60 

22 

82 


Vjerojatnost da student bude žensko i ima ocjenu 3 je 7/82, a 
vjerojatnost student bude muško i ima ocjenu 3 je 15/82. 

Rubna vjerojatnost da bilo koji student ima ocjenu 3 je 22/82, 
a rubna vjerojatnost da student bude muško je 60/82 







Dvodimenzionalne raspodjele 


Za kontinuirane slučajne varijable X,Y definiramo združenu 
funkciju gustoće vjerojatnosti f(x,y), tako da za bilo koji 
dvodimenzionalni intervali c ir x IR vrijedi: 


P((Z,F) E A) = ff f(x,y)dxdy 

A 

Združena funkcija gustoće vjerojatnosti mora zadovoljavati 

uvjete. > o, vx,y II f(x,y)dxdy = 1 


MxR 


Rubne funkcije gustoće vjerojatnosti definiramo kao: 

f xM = f f(x,y)dy frW = / f(?.y)dx 

^ — QO 


Dvodimenzionalne raspodjele 






Primjer: združena funkcija gustoće vjerojatnosti je dana: 


f{x,y) = i 2 +y \x G [0,1],y e [0,1] 


v.0, inače 

Koja je vjerojatnost 
da je X<l/4 i Y<l/4 ? 

Koja je vjerojatnost 
da je X< 1/4 bez obzira 
na Y ? 


f(x,y) 


3.0 



0.50 0.50 


1.00 0.00 
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Uvjetna vjerojatnost 


Neka suXi / kontinuirane slučajne varijable sa združenom 
funkcijom gustoće vjerojatnosti f(x,y) i rubnom funkcijom 
fx(x). Uvjetna funkcija gustoće vjerojatosti za varijablu /, ako 
jeX=x je: 

kvm ‘ MO 

Uz uvjet da je f x (x ) ^ 0 

Analogno pomoću raspodjela vjerojanosti definiramo uvjetnu 
vjerojatnost za diskretne varijable: 

p(x,y ) 




Podsjetnik: nezavisnost događaja 


• Definicija : događaji A i B su nezavisni onda i samo onda 
kad vrijedi: 

P (A | B) = P(A) 

(Vjerojatnost da se dogodi A ne ovisi o tome da li se dogodio B) 

• Iz definicije uvjetne vjerojatnosti slijedi: 

, P(AnB) , N 

PC4|B)= wr =PW 

• Iz toga slijedi: 

P(AHB) = P(A)P(B ) 
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Nezavisne slučajne varijable 

Svaka slučajna varijabla predstavlja ishod nekog događaja, 
stoga možemo poopćiti definiciju nezavisnih događaja na 
nezavisne varijable 

Dvije slučajne varijable će biti nezavisne, ako vrijednost jed 
ne uvjetuje vrijednost druge, tj. ako vrijedi: 

_ p(x,y ) _ 

Pm ~ ~My) ~ Px(x) 

Iz toga slijedi da za nezavisne varijable vrijedi: 

p(x,y) = p x (x)p Y (y ) 



Nezavisne slučajne varijable 


• Općenito, slučajne varijable X i /su nezavisne ako za svaki par 
vrijednosti x i y vrijedi: 


p(x,y) = p x (x)p Y (y ) 


(za diskretne varijable) 


f(x,y) = f x (x)f Y (y) 


(za kontinuirane varijable) 


• Na primjeru raspodjele studenata po ocjeni i po spolu vidimo 
da te varijable nisu nezavisne 
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Nezavisne slučajne varijable 

• Primjer: istovremeno bacanje kocke i novčića 
— slučajna varijabla X opisuje ishod na kocki 
— Slučajna varijabla Y opisuju ishod na novčiću 



• Slijedi da za svaki par vrijednosti x i y vrijedi: 

p(x,y) = p x (x)p Y (y) 

• Dakle varijable su nezavisne! 
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Momenti dvodimenzionalne 

raspodjele 


Pomoćni moment: 




m rs 


X r y s p[x,y') Za diskretne varijable 

xEĐxyED Y 


7Tt rs 


= ff x r y s f(x,y)dxdy 


Za kontinuirane varijable 


MXM 

Centralni moment: 


M rs 


II (x — Mx) r (y — V) Za diskretne varijable 


xED x yED Y 


^TS 


= Jj( x ~ t*x) r (y — iU Y ) S f(x,y)dxdy Za kontinuirane varijable 


RxM 
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Očekivanje 


• Pomoćni momenti m 01 i m 10 : 

m- io = V V x 

x$žDx y£Dy JfG Dx 

m °>=z z 

y£Đy 


yp(x,y ) = V ypy(y) - p Y 


p(x,y) - V xp x (x ) - p x 


• Ovi pomoćni momenti odgovaraju očekivanjima diskretnih 
slučajnih varijabli Xodn. Y 

• Analogno vrijedi i za kontinuirani slučaj 
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Varijanca 


Središnji momenti M 20 i M 02 odgovaraju varijancama 


M 


20 


= y y & - t*x) 2 p(.x>y) 


xED x yED Y 


= ^ (x - Vx) 2 Px(x) = V(X ) 


XED X 


M 


02 


= 7 / (y-py) 2 p0.y) 


xE D x yED Y 


= ^ (y - Py) 2 Py(y) = V(Y) 


yEDy 


Kovarijanca 


• Kovarijanca je veličina koja govori o zavisnosti varijabli X i Y\ 

Cov(X, Y ) = a xy = M n 

• Izražena pomoću pomoćnih momenata ova relacija postaje: 

Cov(X, Y) = m n — m 10 m 01 = E(XY ) - E(X)E(Y) 





tly 


X 


- I + 


• • 

• • 


—r rr _ 

i • • 

+ ! - 


-► x 



Cov(X,Y) > 0 


(i* 

Cov(X,Y) < 0 


Mv 

Cov(X,Y) ~ 0 
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Kovarijanca nezavisnih varijabli 


Po definiciji za kovarijancu vrijedi: 

®xy = E{XY) - E{X)E(Y) 

Prvi pribrojnik možemo pisati kao: 

E(XY) = II xyp(x,y ) 


* y 


Za nezavisne varijable: 

= yy*y Px( x )Py(y ) = V xp x (x) V yp y (y) 


X 


X 


Slijedi: E(XY) = E(X)E(Y ) 

Dobivamo da je kovarijanca nezavisnih varijabli: 


@xy ^ 



Korelacija 

• Kovarijanca ovisi o veličinama koje mjerimo i njihovim jedinicama 
-> nije pogodno za usporedbu povezanosti različitih veličina 

• Uvodimo bezdimenzionalnu veličinu - korelaciju - kao mjeru 
povezanosti, tj. zavisnosti varijabli 

• Koeficijent linearne korelacije varijabli X i Y se definira kao: 


Corr(X,Y ) = p x>Y 


°XY 

o x o Y 


(omjer kovarijance i produkta standardnih devijacija) 


18 




Korelacija 


Primjer 1 : Za 10 studenata uspoređujemo ocjene iz domaće 


zadaće i ocjene na ispitu 


Student 

DZ 

Ispit 

1 

5 

5 

2 

4 

5 

3 

4 

4 

4 

2 

3 

5 

5 

4 

6 

2 

2 

7 

3 

4 

8 

3 

3 

9 

4 

3 

10 

3 

4 


♦ ♦ 

♦ ♦ ♦ 

♦ ♦ ♦ 

♦ 

n-1-1-1-1 

1 2 3 4 5 6 

Ocjena DZ 

C>DZ = 1/02 

^ispit = 0,90 

a Dz-ispu = 0,65 (kovarijanca) 

"^p - 0,7 

Vidimo naznaku linearne ovisnosti 









Korelacija 


• Primjer 2 : Studenti mjere otpor cilindričnog vodiča u ovisnosti 
o njegovom radijusu: 


r/mm 

R/Q 

0,1 

99 

0,2 

26 

0,3 

13 

0,4 

7 

0,5 

3 

0,6 

2,5 

0,7 

2 

0,8 

1,30 

0,9 

1,30 

1 

1,10 


120 
100 
g 80 
O 60 
“ 40 

20 
0 


♦ 


♦ 

♦ 4 

-] -T- T ♦ t ♦ ♦ ♦ ♦ 

0 0,2 0,4 0,6 0,8 1 

r/mm 


1,2 


Koeficijent korelacije p = - 0,68 


Vidimo izraženu obrnuto kvardatnu ovisnost, 


što teorijski i očekujemo: 


R = 


1 l 


onr j 
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Korelacija 


Primjer 3 : promatramo vezu između godine upisa studija i 
trajanja studija 


Godina 

Trajanje /god 

2001 

5,6 

2002 

6,5 

2003 

5,2 

2004 

5,7 

2005 

5,9 

2006 

5,1 

2007 

5,7 

2008 

5,5 

2009 

5,8 

2010 

6,2 


7 -i 

6,5 


(U 

’c 6 

(U 

s_ 

5 5,5 - 

to 
OJ 

i 5 h 


4,5 


♦ 



♦ 


4 -I-,-1-,-,-1-, 

2000 2002 2004 2006 2008 2010 2012 

Godina upisa 


Koeficijent korelacije p = 0,05 
Ne uočava se izražena ovisnost! 







Korelacija 


• Svojstva: 

1. -1 < p < 1 

2. Ako su X i Y nezavisne varijable onda je p=0 -> primjer 3 
(ali obrat ne vrijedi) 

3. Ako je ovisnost varijabli linearna Y=aX+b onda i samo onda je 
p=l ili p=-l (za a>0 odn. a<0) 

4. Može se još pokazati da vrijedi: 

Corr(aX + b,cY + d) = Corr(X, 7) 

• Općenito ako je |p|<0,5 kažemo daje linearna korelacija slaba 
(primjer 3), a ako je |p|>0,5 kažemo da je linearna korelacija 
izražena (primjeri 1 i 2) 

• Iz primjera 2 vidimo da imamo izrazito jaku kvadratnu korelaciju, no 
koeficijent je -0.68 obzirom da je on mjera linearne korelacija 
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Suma slučajnih varijabli 


Definiramo slučajnu varijablu Y=X 1 +X 2 

Očekivanje slučajne varijable Y: 

E(Y) = E(X 1 ) + E(X 2 ) 

Dokaz: 


£■(■^1 + X 2 ) = VV^i + x 2 )p(x 1 ,x 2 ) 

x± x 2 

-II x 1 p(x ll x 2 ) + II x 2 p(x lf x 2 ) 

x± x 2 x ± x 2 

= y x iPxi (^i) +y x 2 p x 2 (x 2 ) 

X± x 2 

= E(X 1 ) + E(X 2 ) 


Suma slučajnih varijabli 

Definiramo slučajnu varijablu Y=X 1 +X 2 

Varijanca slučajne varijable Y: 

V(Y) = V(X 1 ) + V(X 2 )+2 (t x1x2 

Dokaz: 

V(X 1 + X 2 ) = E{[( Xl + x 2 ) - Oi + ,u 2 )] 2 } 

= £{[(x 1 - + (x 2 - ju 2 )] 2 ) 

= E [(Xi - /ii) 2 ] + E[(x 2 - jU 2 ) 2 ] 

+ 2£’[(x 1 - J u 1 )(x 2 - jU 2 )]) 


Linearna kombinacija slučajnih varijabli 

• Za slučajne varijable X v X 2 ,X n možemo definirati novu 
slučajnu varijablu: 

Y = Q ]Xj+Q 2X2i~~‘+Q „X n 

gdje su a v ..a n konstatne. Yse naziva linearnom kombinacijom 
varijabli X. 

• Poopćenje relacije za očekivanje: 

+a 2 Xi H-f a n X n ) = aiE(X\) + a2E{X2')-\ - ha n E(X n ) 


• Vrijedi bez obzira na to da li su varijable međusobno neovisne 
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Linearna kombinacija slučajnih varijabli 


* Poopćenje relacije za varijancu: 



• Posebno, ako su varijable nezavisne, vrijedi: 

V{a 1 X l + a 2 X 2 + --+a n X n ) = afv(X 1 } + ajV(X 2 ) + ■■• + n^V(X„ ) 


* Važan primjer : varijanca razlike nezavisnih varijabli 

nx~Y) = V(X ) + V(Y) = al 

Npr. mjerenje spektra s oduzimanjem pozadine. 
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Linearna kombinacija nezavisnih 
normalnih slučajnih varijabli 


Iz prethodnih razmatranja slijedi specijalni slučaj normalnih 
varijabli: 

Neka su X\, X2,...JC„ nezavisne normalne slučajne varijable. Tada njihova linearna 
kombinacija 

7= aiXi+ a 2 X 2 +...+a n X„ 

n 

također ima normalnu raspodjelu s očekivanjem 1 varijancom 

Z=1 



/=! 


Dokaz preskačemo 




Statistika i osnovna mjerenja 


Teorija uzoraka 

M. Makek 
2016/2017 


Slučajni uzorak 

Iz neke populacije uzmemo uzorak od n - elemenata, koje 
označimo sa x v x n 

U mnogim statističkim problemima se elementi u uzorku 
mogu poistovjetiti sa opaženim vrijednostima n slučajnih 
varijabli X v ... X n 

Definicija : ako su X 1 ,...X n slučajne varijable, one čine slučajni 
uzorak veličine n ako su: 

- Nezavisne 

- Imaju iste raspodjele vjerojatnosti 

(X 1 ,...X n su nezavisne i identično raspodijeljene) 



Slučajni uzorak 

Svaku veličinu koju možemo izračunati iz uzorka nazivamo 
statistikom (npr. prosjek uzorka, standardna devijacija, 
momenti., itd.) 

Obzirom da prije uzimanja uzorka ne možemo znati vrijednost 
statistike ona je također slučajna varijabla 

Primjena : 

— Mjerenja fizikalnih veličina 

— Uzimanje uzorka iz velike populacije 



Prosjek slučajnog uzorka 

1 n 

Je definiran kao: X = — ) X; 

1 = 1 

Prosjek je linearna kombinacija nezavisnih slučajnih varijabli 
X 1 ,... X n te je sam slučajna varijabla 

Očekivanje i varijancu prosjeka možemo odrediti prema 
relacijama za linearnu kombinaciju slučajnih varijabli: 

£(*) = £ (i (*! + ••• + *„)) 

= UE{X ± ) + - + E{X n ))= r ^ i=> 

n n 

vm = v (^(x 1 + -+xS) 

= ^{V{X 1 ') + -+V{X n )) ^ 



V-x = M 






Total slučajnog uzorka 

n 

Je definiran kao: r 0 = %> 

i =1 

Total je linearna kombinacija nezavisnih slučajnih varijabli 
X 1 ,... X n te je sam slučajna varijabla 

Očekivanje i varijancu totala možemo odrediti prema 
relacijama za linearnu kombinaciju slučajnih varijabli: 


£(T 0 ) = i*t 0 = m 

7(T 0 ) = o$ q = no 2 


Raspodjela prosjeka i totala 


Za slučajne varijable s proizvoljnom raspodjelom znamo 
očekivanje i varijancu prosjeka i totala no ne znamo točan 
oblik raspodjele 

Raspodjelu možemo dobiti iz pravila vjerojatnosti (za jedno¬ 
stavne raspodjele) ili pomoću simulacijskog eksperimenta (kad 
je nemoguće provesti analitički račun) 



Raspodjela prosjeka normalne 

raspodjele 

• Specijalni slučaj: normalna raspodjela 

- ako je uzorak sastavljen od n slučajnih varijabli s očekivanjem ju i 
varijancom a 2 , onda i prosjek ima normalnu raspodjelu s očekivanjem 
ju i varijancom d*/n (analogno vrijedi za total) 

- Što je n veći raspodjela prosjeka je uža (što imamo veći broj pokusa 
bolje je definirana srednja vrijednost) 
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Primjer 1: diskretna raspodjela 


Prodavač automobila prodaje 50% automobila niže klase po cijeni 65.000 
kuna, 35% automobila srednje klase po 100.000 kuna i 15% automobila 
vise klase po 135.000 kuna. 

Definirajmo slučajnu varijablu 

X=prihod od prodaje jednog automobila (u tisućama kuna) 



0 . 8 -, 


fi =87,75 

<7 2 =646,19 

<t=25,42 


65 0,5 
100 0,35 
135 0,15 



0 . 6 - 


30 


65 


100 


135 


X 
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Primjer 1: diskretna raspodjela 

Određenog dana najavila su se dva kupca. Neka su slučajne varijable: 
X x =pnhod od prvog kupca 
X 2 =prihod od drugog kupca 


Tablica mogućih prihoda: 


Xi 

X 2 

lXXl.X2> 

Xi+X2 


65 

65 

0.25 

130 

65 

65 

100 

0.175 

165 

82.5 

65 

135 

0.075 

200 

10) 

100 

65 

0.175 

165 

82.5 

100 

100 

0.1225 

200 

100 

100 

135 

0.0525 

235 

117.5 

135 

65 

0.075 

200 

100 

135 

m 

0.0525 

235 

117.5 

135 

135 

0.0225 

270 

135 


Slučajna varijabla X 
raspodijeljena je ovako: 


X 

p(x) 

65 

0,25 

82,5 

0.35 

100 

0,2725 

117,5 

0.105 

135 

0,0225 



X 
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Primjer 1: diskretna raspodjela 





X 
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Primjer 1: diskretna raspodjela 



X 


lg 

CL 


D,4 -i 


0,3- 


D,2 - 


□ . 1 - 


0.0 


025 - 


= 87,75 
cr 2 = 323,09 



3D 65 1DD 135 


020 - 



0.10 - 


0.05 - 


0.00 


X 



X 


n = 1 


n = 2 


n = 4 


ii 












Primjer 2: kontinuirana raspodjela 

Na putu do posla čekam autobus koji vozi svakih 10 minuta, a zatim tramvaj 
koji također vozi svakih 10 minuta. 

Neka su slučajne varijable: 

X 1 = vrijeme čekanja autobusa 
X 2 = vrijeme čekanja tramvaja 

1/10 , 0 < x -> <10 

f(x 2 ) = *1 

' 1 0 , X 2 <0 ili X 2 >i0 


/(u) 


1/10 , 0 < ^ <10 
0 , x\ < 0 ili > 10 


_ 0 / 1 - 

XT 
sk r 


0,1 


X 


- 

0 


r -T-,- T -T- T -T-,- ---1-1 .-- 0,0 -T-,-^ 

10 0 

T 0 = ukupno vrijeme čekanja =X l +X 2 


10 
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Primjer 2: kontinuirana raspodjela 



+ 




t 


m = 


100 
20- t 

100 
0 , 




0 < t < 10 

,10 < t < 20 
inače 
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Uočiti 

Iz gornjih primjera vidi se da čak i kada populacija 
nije normalno raspodijeljena, raspodjela prosjeka 
totala teži normalnoj raspodjeli kada je uzorak 
dovoljno velik 

-> Centralni granični teorem 



Centralni granični teorem 

• Neka je X 1/ ... X n slučajni uzorak proizvoljne raspodjele 
očekivanjem //i standardnom devijacijom a. Ako je n dovoljno 

velik X ima približno normalnu raspodjelu s očekivanjem Px = t* 

2 2 

i standardnom devijacijom % = ° i n 

7 0 ima približno normalnu raspodjelu s očekivanjem Mr 0 — 

i standardnom devijacijom oi = no 2 

• U praksi se teorem može primijeniti za n>30 
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Centralni granični teorem 


• Primjer: 



• Simulacija: 

http://onlinestatbook.com/stat_sim/sampling_dist/ 
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Statistika i osnovna mjerenja 


Procjenitelji i pogreške 

M. Makek 
2016/2017 


Procjena parametara populacije 


Uzmimo da populacija ima karakteristične parametre poput /u, 
a, o 2 , momenti, itd. 

Te parametre je nemoguće odrediti na cijeloj populaciji, stoga 
se uzima uzorak i na temelju njega se procjenjuju parametri 
populacije 

Uzmimo slučajnu varijablu © kao procjenitelj parametra 0 
©je nepristrani procjenitelj parametra 0 ako vrijedi: 

g(0) = 9 

tj. mu je očekivanje stvarno jednako tom parametru 



Procjena parametara populacije 


Ako je X v X n slučajni uzorak populacije, onda je slučajna 
varijablažnepristrani procjenitelj očekivanja populacije jj, 


Dokaz : pok azali smo da je o čekivanje^ jednako očekivanju 


populacije E(X ) = = fj. 


Ako je X 1 ,X n slučajni uzorak populacije, onda je slučajna 
varijabla S 2 nepristrani procjenitelj varijance populacije <y 2 



Dokaz slijedi 






Procjena parametara populacije 


Pokazali smo da za bilo koju slučajnu varijablu V vrijedi: 

V(Y)=E(Y 2 )- [J E(Y )] 2 

Stoga za slučajnu varijablu X,- možemo pisati: 

E(X?) = [E(X t )] 2 + V(X t ) 

Raspišemo izraz za S 2 : 

52 -;£*■)' 


2i 




Procjena parametara populacije 


B(s2) =;rhE B «> - ME *<)1 j 



-frE*M'G> 

= ——r|V (ju 2 + * 2 ) - ~{V[T 0 ] + [E(T 0 )] 2 }- 

n — 1 l/—t n J 

E (ju 2 + <r 2 )-{ncr 2 + [n^] 2 } = <t 

n ; 


Što je i trebalo pokazati 








Intervali pouzdanosti 


Procjenom parametara dobivamo njihove iznose, no ne znamo 
koliko su ti iznosi pouzdani! 

Npr. ako izračunamo prosjek uzorka on neće biti identičan 
očekivanju populacije, no zanima nas koliko je njihovo 
odstupanje 

Prema središnjem graničnom teoremu, srednja vrijednost 
uzorka ima normalnu raspodjelu X~N(n, a 2 /ri) gdje su // i cr 
očekivanje i varijanca populacije 


Pouzdanost procjene prosjeka 

populacije 


Obzirom da srednja vrijednost uzorka ima normalnu 
raspodjelu, vrijedi: 

p 


, x - u , 
-i<— h=<i 


\ 


<j 




J 


= 68 % 


X- LI 
-1.96 < —t-f= < 1,96 
<x/ 




= 95% 




- 2.575 < , ^ < 2.575 




nj-Ju 


= 99% 




Dakle možemo reći da se prava vrijednost očekivanja 
populacije nalazu i intervalu: 

s vjerojatnošću od 68% 













Pouzdanost procjene prosjeka 

populacije 


u prethodnoj relaciji smo pretpostavili da je varijanca 
populacije a poznata. 


Ako nije poznata onda je procjenjujemo pomoću slučajne 
varijable: 


5 = 


- x) 

n— 1 


Tada interval od 68% pouzdanosti postaje: 





Pogreške mjerenja 


U fizici se rezultati mjerenja pišu kao srednj a vrijed nost nezavisnih 


mjerenja u intervalu od 68% pouzdanosti: x ± M 


Gdje je 


nazivamo i 




M = 


n(n — 1 ) 


nepouzdanost prosjeka, koju 


standardnom greškom (sad zn amo od kuda dolazi ova form ula) 

Pri tome je rasipanje odn. 

Standardna devijacija uzorka: m = s = 



Obzirom da je rasipanje dolazi zbog nemogućnosti kontrole slučajnih 
pogrešaka, kažemo da je m mjera preciznosti aparature 

m se ne mijenja značajno s povećanjem broja mjerenja , no 
M=m/Vn se smanjuje te kažemo da je rezultat pouzdaniji ako 
napravimo više mjerenja 









Mjerenja različitih statističkih težina 

• Mjerimo neku fizikalnu veličinu u nekoliko različitih mjerenja s 
različitim pouzdanostima zanima nas procjena prave 
vrijednosti te veličine na temelju tih mjerenja 

• Izvedeno je k nizova mjerenja te fizikalne veličine te je za svaki 
niz dobivena srednja vrijednost i standardna pogreška: 

x t = x t ± M ir i = 1... k 

• Tražimo najvjerojatniju vrijednost x„ mjerene fizikalne veličine 
na temelju tih mjerenja 
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Opća srednja vrijednost 


Funkcija gustoće vjerojatnosti je na temelju k serija mjerenja 

1 


f i > * * ■ j “^k ) 


M 1 ... M/^VŽžr) 


_ yfe fe ) 
1=1 2 M? 


Da bi našli za koji x p je/maksimalan logaritmiramo i tražimo 


minimum sume: 


k 


1 


i = 1 


2M} 


— min. 






Opća srednja vrijednost 


Iz uvjeta minimuma slijedi: 

d \ 1 

dx Ti i—L 

y i =1 


k 


Zj t =l 


2Mf 


x, 


-I 


( x i x p ) 


M; 


= 0 


i=l 


M; 


Xp = 


“i=i 


Pri tome definiramo statističke 
težine pojedine serije kao: _L 


M; 


Mf 


Wi - 


Slijedi opća srednja vrijednost: 


y« 

1 


O 



M, 2 


12 











Nepouzdanost opće srednje vrijednosti 


Da bi odredili pouzdanost sjetimo se da je procjenitelj prave 
vrijednosti veličine, linearna kombinacija: 


k 


x, 


= ^ w i X i 


i=l 


Prema relacij i za varijancu l inearne kombinacije dobivamo: 


k 

k 

m 2 =v(x p )= y 

wfV{X i ) = V wf M 2 

i -1 

i =i 


Primjetimo još da se statističke težine onose kao: 

11 1 

w 1 :w 2 : w n — 


I da vrijedi: 


Ml'Ml . Ml 


k 


'Z w ‘= 


i =1 
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Propagacija pogreške 

• Želimo odrediti neku fizikalnu veličinu preko njene funkcijske 
ovisnosti u neposredno mjerenim veličinama 

• Primjeri: 

a) Mjerimo masu i volumen predmeta i želimo odrediti njegovu gustoću 

b) Mjerimo put i vrijeme koje auto prijeđe i želimo odrediti njegovu brzinu 

• Pitanje je kako pogreška posredno mjerene veličine ovisi o 
pogreškama neposredno mjerenih veličina? 

• Pri sljedećim razmatranjima pretpostavljamo da su pogreške 
mjerenja slučajne (normalno raspodijeljene) i malene tj. da 
vrijedi Mj«jUj 

• Razmotrimo prvo dva specifična primjera koja ćemo poopćiti 
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Propagacija pogreške 

1. Posredna veličina /je linearna kombinacija izravno mjerenih 
veličina X i Y: f(X,Y)=aX+bY 

-> Za slučaj linearne kombinacije nezavisnih varijabli znamo 
izračunati očekivanje i varijancu, dakle: 

E[f(X, r)] = aEO o + bE(Y) M f = a 2 Mf + b 2 M$ 

2. Posredna veličina /je nelinearna funkcija veličine X: f(X) 

-> za očekivanje vrijedi (uz pretpostavku malene pogreške): 

E[f(X)] ^ f{pL x ) = f(X) 

-> za varijancu vrijedi: 

Mf = V[f(X )] = E [(/(X) - /( fi x )) 2 
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Propagacija pogreške 


dalje razvijemo f(X) u red: f (X) = /(&?) + (X -1*x ) 


df 


dX 


+ 


slijedi da je: 


E [(/(*) - - E 


(X - ii x ) 


df 


dX 




l*x 


df 


dX 




'v-xj df 

odnosno: f(X) - f(p x ) - (X - (i x ) — 

đX 


Dakle za male pogreške vrijedi: M f = 


f*x 


df 


dX 



M, 


16 
















Propagacija pogreške 

3. Općeniti slučaj - posredna veličina /je funkcija n izravno 
mjerenih neposrednih nezavisnih veličina X v X n 

-> pokazuje se da tada vrijedi: 


/(*!, . 

■ ■ / ^TI ) f(%l> t 


II 

li®, ,)'* 

t=l v x 1 ) ...x n / 


To su relacije koje smo uveli na početku, a sada je objašnjeno 
odakle one dolaze 
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Propagacija pogrešaka s korelacijom 

• Prethodne relacije za propagaciju pogreške vrijede kada su 
slučajne varijable X i Y nezavisne 

• Sada promatramo općenit slučaja kada dopuštamo da 
mjerenja X i Y nisu međusobno nezavisna 

• Uzmimo da smo izvršili n mjerenja X i Y i kao rezultat dobili 
{x,y,) t gdje i=l...n. Na temelju dobivenih vrijednosti 
izračunamo n vrijednosti veličine q i =q(x i ,y i ) 

• Zanima nas srednja vrijednost i pogreška veličine q 
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Propagacija pogrešaka s korelacijom 


Pretpostavljamo da su pogreške opservabli X i Y relativno male 
tj. da su odstupanja x t — x i y; — y mala tako da pri razvoju 
u red možemo uzeti samo linearne članove: 




dq 

q(*,y) + Oi-*) 


dx 


dq 

+ C yi - y) t- 

x,y oy 


x t y 


Prema definiciji srednja vrijednost je: 


, n n r 

i = l 1=1 


dq 

q(x,y) + (x t - x) — 

s _ 


T 

0 


, , -,dq 

_ 


x,yi 


Slijedi: q = q (x, y) 


o 
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Propagacija pogrešaka s korelacijom 


Prema definiciji standardnu devijaciju možemo dobiti prema: 

n n r- T 


= nZ/’ 1 ‘ ^ = sZ 


i=1 
n 


i= 1 




x f y 


,, -,dq | 
+ O',-) 05 d 




- (D’;Ž* - +©’£<* - ** -» 

£=1 i=l y 1=1 



• gdje je (j kova rija nca 

• Gornja relacija daje standardnu devijaciju veličine q bez obzira 
na to jesu li veličine X i Y povezane. 

• Za nezavisne varijable X i Y kovarijanca je 0, a relacija za 
standardnu devijaciju se svodi na prva dva člana 
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Propagacija pogrešaka s korelacijom 


• Ako kovarijanca ne iščezava kažemo da su pogreške opservabli 
Xi Vkorelirane, a standardna devijacija cr q je očito različita 
nego u slučaju neovisnih, slučajnih pogrešaka. 






Može se pokazati da relacija za <j q daje opću gornju granicu 
standardne devijacije veličine q. 


Da bi to dokazali iskoristit ćemo Cauchy-Schwarzovu 
nejednakost: 


®xy — ®x®y 


Uvrštavanjem te nejednakosti u izraz za standardnu devijaciju: 



< 



2 

+ 



(Ty + 2 


dq dq 
dx dy 


®x ®y 

x f y 
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Propagacija pogrešaka s korelacijom 


Dalje slijedi: 


r 2 

q 


< 


dq 

(7 + 

dq 

(Ti, 



dx 

X 1 

dy 



.,2 


Čime dobivamo gornju granicu za standardnu devijaciju bez 
obzira na korelaciju varijabli X i Y 


a n < 

dq 

O v + 

dq 

<7^, 

q — 

dx 

X 1 

dy 

y 
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Statistika i osnovna mjerenja 


Testiranje hipoteza 

M. Makek 
2016/2017 


Postavljanje hipoteza 


• Postavljamo neku teoriju i želimo provjeriti da li je ona istinita 


• Primjeri : 

1. Kocka (za bacanje) nije poštena 

2. Novi lijek je bolji od starog 

3. Teorijski model dobro opisuje izmjerene podatke 

4. Optuženi nije kriv 


• Da bi se izvršila provjera postavljamo dvije hipoteze: 


a) H 0 - nul-hipoteza 

b) Hj- alternativna hipoteza 


-> Samo jedna od njih je točna! 


• H 0 se smatra ispravnom dok se ne dokaže suprotno H 0 i 
se ne tretiraju ravnopravno, već se konzervativnija postavlja 
kao istinita dok se ne dokaže suprotno. 

Ako ne odbacimo H& to ne znači da je ona ispravna, nego samo da nemamo dovoljno 
dokaza da ju odbacimo. 




Postavljanje problema 


1. Prvo određujemo: 

- A - područje prihvaćanja hipoteze H 0 

- B - Područje odbacivanja hipoteze H 0 (kritično područje) 


Uzimamo uzorak iz populacije koju istražujemo 

n An Ar I rA a Ari 1 1 i 1/ 1 i • 

L/UI IUjI i i i vj uuiu ivu . 

ODLUKA 


Odbaci H 0 

Prihvati H 0 

ISTINA H ° 

Hl 

Pogreška 1. vrste 

Ispravan zaključak 

Ispravan zaključak 

Pogreška 2. vrste 


4. Određujemo vjerojatnost pogreške: 

- a= P(B | H 0 ), vjerojatnost da odbacimo H 0 kada je istinita 

- (3 = P(A | H!), vjerojatnost da prihvatimo H 0 kada je neistinita 

-> Smatra se da je pogreška 1. vrste teža pogreška (npr. krivnja 
optuženika) 



Signifikantnost testa 


Provjera (test) ima signifikantnost a, ako je vjerojatnost za 
pogrešku prve vrste < a. (Test razine a) 

Konvencionalno se za razinu signifikantnosti uzima: 

— a=0,05 (signifikantan) 

— a=0,01 (vrlo signifikantan) 

Snaga testa je vjerojatnost da test uputi na ispravnu odluku, 
tj. da odbacimo H 0 kad je uistinu neistinita: P=l-|3 



Primjer - bacanje novčića 

Hipoteze: 

— H 0 = novčić je pošten 

— H 1 = novčić nije pošten 

X = broj 'pisama' u 6 bacanja -> ako je novčić pošten onda je 
vjerojatnost P(X | H 0 ) ~ Bin(6; 0.5) 

Mogući su sljedeći ishodi: 



Odaberemo: A za X=l,2,3,4,5 i B za X=0,6 

Vjerojatnost pogreške 1. vrste a=P(B|H 0 )=0,0156 x 2 

Dakle ako u 6 bacanja kocke dobijemo 0 ili 6 puta pismo i na temelju 
toga zaključimo da je novčić nepošten, vjerojatnost da smo donjeli 
pogrešnu odluku je 3,12% 









Prilagodba krivulje mjerenim 

podacima 

Pretpostavimo da smo izmjerili neku raspodjelu događaja (raspodjelu 
frekvencija) 

Prema obliku raspodjele i teoriji 
(ako je poznata) možemo mjerenim 
podacima prilagoditi krivulju regre¬ 
sijskim metodama (pravac, polinom, 

Gaussijan, Poissonova raspodijela, 
binomna raspodjela, itd.) 

X 

Obzirom da svako mjerenje ima ograničenu preciznost javit će se 
odstupanja od krivulje 

Pitamo se da li su ta odstupanja slučajne prirode i da li pretpostavljena 
krivulja dobro opisuje izmjerene podatke 





Dobrota prilagodbe 


Da bi ispitali kvalitetu (dobrotu) prilagodbe izmjerenim 
podacima izvodimo statistički test 

Jedan od testova koji se najčešće koriste je % 2 -test 

Podsjetimo se T— funkcije, općenite T- raspodjele i specifično 
X 2 raspodjele 


Gama funkcija 


Definiramo T-funkciju: 

• Svojstva : 

1. Va > 1 -> r(a) = (a — l)r(a — 1) (rekurzija) 

2. Vn £ M -> T(n) = (n — 1)! (poopćenje faktorijela!) 


co 

r(a) — f x a ~ 1 e~ x dx, 

a > 0 

J 

0 




30 n 









Općenita T-raspodjela 

• Funkcija gustoće vjerojatnosti: 

{ X a-l e -x/p 

p«r(a) 

0 , 

gdje su a, (3 >0 

• Očekivanje: E(X)=aj3 
• Varijanca: V(X)=a/3 2 


x > 0 
x < 0 



% 2 -raspodjela 


Specijalni slučaj općenite T raspodjele za (3=2 i a=v/2 
Funkcija gustoće vjerojatnosti: 


r x v/2-l e -x/2 


x > 0 


f(x;v)= < 2 v / 2 r(v/2) ' 

V 0, x < 0 

gdje je v prirodan broj i naziva se "broj stupnjeva slobode" 

• Veza s normalnom raspodjelom : 

Tm. Ako je u varijabla koja ima normalnu raspodjelu s očekivanjem 
standardnom devijacijom ju i o, onda je varijabla 


x = 


/ \ 2 

u-JU 

V cr J 


raspodijeljena po % 2 -raspodjeli. Dokaz preskačemo. 







X 2 -test 


1. Uzmimo da eksperiment daje raspodjelu frekvencija//, koje 
grupiramo u n razreda (i=l...n). 

2. Izmjerenim podacima pomoću regresije prilagodimo krivulju, koja 
daje teorijske frekvencije za pojedini razred f ti 



3. Definiramo veličinu: 


4. Paersonov teorem : veličina tfje približno raspodijeljena prema x*~ 
raspodjeli sa v stupnjeva slobode koji ovise o broju opažanja (n) i 
broju ograničenja (k): v= n-k 

5. Pri tome diskretnu raspodjelu aproksimiramo kontinuiranom. 
Aproksimacija ne vrijedi ako je/ t/ < 5 -> u tom slučaju treba 
grupirati razrede tako da je frekvencija pojedinog razreda uvijek 
veća od 4 


ii 




Signifikantnost x 2 -testa 


• Ako su razlike |f r f ti | relativno male, možemo zaključiti da imaju slučajan 
karakter 


• Ako su razlike |f r f ti | relativno velike (signifikantne), onda će i % 2 biti velik 
(signifikantan), tj. možemo zaključiti da je prevelik da bi bio slučajan. Tada 
je razumno odbaciti hipotezu. 

• Vjerojatnost da varijabla X poprimi vrijednost veću od neke zadane Xv( a )> 
jednaka je površini repa ispod krivulje i dana je integralom: 


P(X > Xv( a )) = jffa v) dx 

Xl(a) 

• Po dogovoru se uzima da je % 2 signifikantan 
ako padne u kritičnih 5% raspodjele i 

u tom slučaju odbacujemo hipotezu 

• U prosječno 5% slučajeva će odbacivanje 
biti pogrešno jer je u toliko slučajeva 
moguće da je x 2 signifikantan, a hipoteza 
je istinita. 
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Primjer 1 


Neki uređaj daje određeni postotak neispravnih proizvoda. 
Proizvodi se pakiraju u kutije od 20 komada. Za probu se 
otvara 50 kutija i određuje se broj neispravnih proizvoda X 


Opaženo je sljedeće: 



srednji x=l,8 

Koja je vjerojatnost da se u kutiji nađe X defektnih komada? 

— teorijski bi trebala biti dana binomnom raspodjelom: za x = np i 
n=20, dobivamo da je p=0.09 


, , / 20 \ 
h*) = M 


x 20 —x 

pq 


- Teorijske frekvencije: 
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Primjer 1 


Nul-hipoteza: binomna raspodjela uz pretpostavku od 9% 
defektnih proizvoda dobro opisuje podatke 

Želimo potvrditi/odbaciti H 0 uz signifikantnost 0.05 

Broj stupnjeva slobode: 

— Broj razreda n=5 

— Broj ograničenja: 2 (određen p i N=50) 

— v=5-2=3 

Izračunamo % 2 = 3,08 

Iz tablica vidimo da je za 3 stupnja slobode i 5% 
signifikantnosti X 2 o,os ; 3 = 7/82 

Obzirom da dobivena vrijednost ne spada u kritično područj 
tj. x 2 < x 2 0 , 05 ;b zadržavamo nul hipotezu. 



Primjer 2 


Dugogodišnje statistike pokazuju daje visina studenata na nekom 
sveučilištu normalno raspodijeljena s očekivanjem // = 173 cm i 
standardnom devijacijom a=l cm. 

Iz jedne generacije studenata 11 a tom sveučilištu izdvojeno je 100 
studenata, mjerene su njihove visine i svrstane u razrede. 

Opažene su frekvencije: 


X 

154-160 

160-166 

166-172 

172-178 

178-184 

184-190 

Total 

f(x) 

5 

17 

38 

25 

9 

6 

100 


Provjeri s 5% signifikantnosti je li raspodjela normalna s s 
očekivanjem // = 173 cm i standardnom devijacijom a= 1 cm ! 


15 










Primjer 2 

H 0 : Raspodjela je Gaussova s očekivanjem jn = 173 cm i 
standardnom devijacijom a = 1 cm. 



visina 


Imamo jedno ograničenje: lf i =If ti =N= 100 
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Primjer 2 


Izmjerene i teorijske frekvencije su dane tablicom: 


razred i 

f, 

ft, 

(frfnf/ftt 

154-160 

5 

3 


160-166 

17 

13 

22 5 

166-172 

QO 

28 

3,57 

172-178 

25 

32 

1,53 

178-184 

9 

18 

4,5 

184-190 

6 

6 

0 

zbroj 

100 

100 

i 1,85 


Tablica ima 5 razreda te imamo jedno ograničenje N=100 
->broj stupnjeva slobode je 4 












Primjer 2 


Za 4 stupnja slobode i 5% signifikantnosti, kritična vrijednost 
(iz tablica) je 0 

Zo. 05,4 =9.49 



Izračunata vrijednost /op = 11-85 pada u kritično područje. 
Stoga hipotezu H 0 odbacujemo. 
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Primjer 2 


• Nova hipoteza: 

H 0 : Raspodjela je Ganssova s očekivanjem p = -^uzorka 
i standardnom devijacijom <J = cr uzor k a . 

Izračunamo: 

- i n i 

x uzorka = ^ fi x i = 1 ^ 1 cm ^uzorka = J jj S fi ( x ? — x uzorka ) = ^,1 cm 



visina 
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Primjer 2 


• Izmjerene i teorijske frekvencije su dane tablicom: 


razred /' 

f 

ft, 

(fr.fi) 2 (fti 

154-160 

5 

6 

0,167 

160-166 

17 

18 

0,056 

166-172 

38 

32 

1,125 

172-178 

25 

28 

0,321 

178-184 

9 

13 

0,063 

184-190 

6 

3 


zbroj 

100 

100 

1,732 


Dakle, imamo pet razreda i tri ograničenja: N= 100 

// — * uzor k a 


® ^uzorka 


pa je v = 5-3 = 2 
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Primjer 2 


Za 2 stupnja slobode i 5% signifikantnosti, kritična vrijednost j 

Zo.05,2 = 



7 

Izračunata vrijednost Xop = b73 ne pada u kritično područje. 
Stoga hipotezu H 0 zadržavamo. 





Dodatak: % 2 -tablica 


df 

X 2 93S 

X 2 D90 

X%75 

X%50 


X 2 100 

X 2 05O 

X 2 025 

\Joio 


1 

0.000 

0.000 

0.001 

0.004 

0.016 

2.706 

3.841 

5.024 

6.635 

7.879 

2 

0.010 

0.020 

0.051 

0.103 

0.211 

4.605 

5.991 

7.378 

9.210 

10.597 

3 

0.072 

0.115 

0.216 

0.352 

0.584 

6.251 

7.815 

9.348 

11.345 

12.838 

4 

0.207 

0.297 

0.484 

0.711 

1.064 

7.779 

9.488 

11.143 

13.277 

14.860 

5 

0.412 

0.554 

0.831 

1.145 

1.610 

9.236 

11.070 

12.833 

15.086 

16.750 

6 

0.676 

0.872 

1.237 

1.635 

2.204 

10.645 

12.592 

14.449 

16.812 

18.548 

7 

0.089 

1.239 

1.690 

2.167 

2.833 

12.017 

14.067 

16.013 

18.475 

20.278 

8 

1.344 

1.646 

2.180 

2.733 

3.490 

13.362 

15.507 

17.535 

20.090 

21.955 

9 

1.735 

2.088 

2.700 

3.325 

4.168 

14.684 

16.919 

19.023 

21.666 

23.589 

10 

2.156 

2.558 

3.247 

3.940 

4.865 

15.987 

18.307 

20.483 

23.209 

25.188 

11 

2.603 

3.053 

3.816 

4.575 

5.578 

17.275 

19.675 

21.920 

24.725 

26.757 

12 

3.074 

3.571 

4.404 

5.226 

6.304 

18.549 

21.026 

23.337 

26.217 

28.300 

13 

3.565 

4.107 

5.009 

5.892 

7.042 

19.812 

22.362 

24.736 

27.688 

29.819 

14 

4.075 

4.660 

5.629 

6.571 

7.790 

21.064 

23.685 

26.119 

29.141 

31.319 

15 

4.601 

5.229 

6.262 

7.261 

8.547 

22.307 

24.996 

27.488 

30.578 

32.801 

16 

5.142 

5.812 

6.908 

7.962 

9.312 

23.542 

26.296 

28.845 

32.000 

34.267 

17 

5.607 

6.408 

7.564 

8.672 

10.085 

24.769 

27.587 

30.191 

33.409 

35.718 

18 

6.265 

7.015 

8.231 

9.390 

10.865 

25.989 

28.869 

31.526 

34.805 

37.156 

19 

6.844 

7.633 

8.907 

10.117 

11.651 

27.204 

30.144 

32.852 

36.191 

38.582 

20 

7.434 

8.260 

9.591 

10.851 

12.443 

28.412 

31.410 

34.170 

37.566 

39.997 


22 


















